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We present a new (semi-) analytic model for feedback in galaxy formation. The 
interstellar medium (hereafter ISM) is modeled as a two-phase medium in pressure 
equilibrium, where the cold phase is fragmented into clouds with a given mass spec- 
trum. Cold gas infalls from an external halo. Large clouds are continually formed 
by coagulation and destroyed by gravitational collapse. Stars form in the collapsing 
^ . clouds; the remnants of exploding type II supernovae (hereafter SNe) percolate into 

ly-^ ' a single super-bubble (hereafter SB) that sweeps the ISM, heating the hot phase (if 

T-ij- I the SB is adiabatic) or cooling it (in the snowplow stage, when the interior gas of the 

CnI ' SB has cooled). Different feedback regimes are obtained whenever SBs are stopped 

' either in the adiabatic or in the snowplow stage, either by pressure confinement or by 

, blow-out. 
JnI ■ The resulting feedback regimes occur in well-defined regions of the space defined 

by vertical scale-length and surface density of the structure. In the adiabatic blow-out 
(~| . regime the efficiency of SNe in heating the ISM is rather low (~5 per cent, with ~80 

i' per cent of the energy budget injected into the external halo), and the outcoming ISM 

' I is self- regulated to a state that, in conditions typical of our galaxy, is similar to that 

^ . found in the Milky Way. Feedback is most efficient in the adiabatic confinement regime, 

where star-formation is hampered by the very high thermal pressure and the resulting 
^ , inefficient coagulation. In some significant regions of the parameter space confinement 

takes place in the snowplow stage; in this case the hot phase has a lower temperature 
and star formation is quicker. In some critical cases, found at different densities in 
several regions of the parameter space, the hot phase is strongly depleted and the cold 
;J] ' phase percolates the whole volume, giving rise to a burst of star formation. 

I While the hot phase is allowed to leak out of the star-forming region, and may 

give rise to a tenuous wind that escapes the potential well of a small galactic halo, 
strong galactic winds are predicted to happen only in critical cases or in the snowplow 
confinement regime whenever the SBs are able to percolate the volume. 

This model provides a starting point for constructing a realistic grid of feedback 
solutions to be used in galaxy formation codes, either semi-analytic or numeric. The 
predictive power of this model extends to many properties of the ISM, so that most 
parameters can be constrained by reproducing the main properties of the Milky Way. 

Key virords: galaxies: formation - galaxies: ISM - ISM: bubbles - ISM: kinematics 
and dynamics 



1 INTRODUCTION 

Galaxy formation is an open problem. This is due to the 
complexity of the feedback processes that arise from the 
energetic activity of massive or dying stars, taking place 
through winds, ionizing photons and SN explosions (not to 
mention AGN). These feedback processes involve a largo 
range of scales and masses, from the sub-pc scale of star 
formation to the >10 kpc scale of galactic winds, and from 
1 to 10^^ Mq or more. 



It is useful at this stage to identify ranges of scales 
in which different processes are dominant. On >1 kpc spa- 
tial and >10^ Mq mass scales the dominant processes such 
as shock heating of gas, radiative cooling, disc formation, 
galaxy merging and tidal or ram-pressure stripping are 
closely related to the dark-matter halo hosting the galaxy 
and to its hierarchical assembly. On scales ranging from ~1 
pc to ~1 kpc, or from ~1000 to ~10^ Mq, cool gas reaches 
suitable conditions for collapse and star formation, and the 
energy input from massive stars (through winds, UV pho- 
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tons and SNe) acts in shaping and sustaining the muhi- 
phase structure of the ISM. At smaller scales star forma- 
tion takes place; it is most likely driven and self-limited by 
magucto-hydro-dynamical (MHD) turbulence. This division 
is obviously meant to be only a rough approximation of re- 
ality. 

Numerical simulations of whole galaxies are still limited 
to space and mass resolutions not much smaller than ~1 kpc 
and 10^ Mq respectively (see, e.g., Weinberg, Hernquist & 
Katz 2002; Steinmetz & Navarro 2002; Mathis et al. 2002; 
Lia, Portinari & Carraro 2002; Recchi et al. 2002; Peaxce et 
al. 2001; Toft et al. 2002; Springel & Hernquist 2003; Tor- 
natore et al. 2003; Governato et al. 2004). They can address 
effectively the processes dominant in the large-scale range 
identified above, but the feedback processes acting on in- 
termediate and small scales are "sub-grid" physics and are 
treated with simple heuristic models that require the intro- 
duction of free parameters. 

Current models of semi-analytic galaxy formation treat 
feedback at a similar, phenomenological level (see, e.g.. Cole 
et al. 2000; Somerville, Primack & Faber 2001; Diaferio et al. 
2001; Poh et al. 2001; Hatton et al. 2003); they typically con- 
nect the efficiency of feedback to the circular velocity of the 
dark matter halo, with the aid of free parameters. Models of 
galaxy formation that include a more detailed description of 
feedback have been presented, e.g., by Silk (1997, 2001), Fer- 
rara & Tolstoy (2000), Efstathiou (2000), Tan (2000), Lin & 
Murray (2000), Hirashita, Burkort & Takcuchi (2001), Fcr- 
reras, Scannapieco & Silk (2002) or Shu, Mo & Mao (2003). 
In this framework (semi-)analytic work can give a very use- 
ful contribution in selecting the physical processes that are 
most likely to contribute to feedback. 

The focus of this paper is on modeling the intermediate 
range of scales defined above, where the physics of the ISM 
is in act. The standard picture of the ISM is that of a multi- 
phase medium in rough pressure equilibrium; the reference 
model is that of McKcc & Ostrikor (1977), who considered a 
medium composed by cold, spherical clouds with tempera- 
ture and density Tc^lOO K and ric^lO cm~^, kept confined 
by a hot phase with Th~10® K and nh~10~^ cm~^. A warm 
phase of r„~10'* K and nw-^lO"^ cm~^ was produced at 
the interface. This vision is partially confirmed by multi- 
wavelength observations (see, e.g., Heiles 2001), although 
reality appears more complex, suggesting the presence of at 
least 5 different phases. 

This picture is challenged by the results of many sim- 
ulation programs, aimed to the numerical modeling of the 
ISM (see, e.g., Mac Low et al. 1998; Ostriker, Gammie & 
Stone 1999; Avila-Rooso & Vazquez-Somadeni 2001; Kritsuk 
& Norman 2002; see Mac Low 2003 and Vazquez- Semadeni 
2002 for reviews). In this context the ISM is dominated by 
compressible, supersonic, MHD turbulence. These groups 
are still struggling to tame the full complexity of the prob- 
lem, so that these simulations are not directly aimed to 
or easily usable by modeling of galaxy formation. For our 
purposes it is worth mentioning some results. The distribu- 
tions of temperature and density of the simulated gas parti- 
cles show a wide range of values without any strong multi- 
modality, but some broad peaks arc anyway present. The 
distribution of pressure shows a much more limited range 
of values. Structures defined as overdensities are not static 
clouds but transient features of an overall fractal distribu- 



tion (which is consistent with observations, see Chappell & 
Scale 2001) that do not last more than a sound crossing time, 
unless they are gravitationally bound. Thus, the "classical" 
picture of the ISM is not validated, but a model with mul- 
tiple phases in rough pressure equilibrium can still be used, 
though with care, as a useful first-order approximation, able 
to catch some significant elements of the dynamics of the 
ISM. 

The motivation for the present work is to investigate 
the kind of physical processes that arise in galaxy forma- 
tion, in order to provide a grid of solutions for the behaviour 
of feedback in a wide range of realistic cases, to be used 
in simulations or semi-analytic models of galaxy formation. 
We restrict to a two-phase medium in pressure equilibrium, 
composed by cold clouds embedded in a diffuse hot phase. 
The dynamics of the ISM is at present Eissumed to depend 
only on its "local" properties, leaving thus out "large scale" 
events like differential rotation, spiral arms, mergers, galac- 
tic winds and so on. These events will be introduced once 
the global characteristics of the galaxy are specified. 

This paper is the first of a series aimed to modeling 
feedback in galaxy formation. It presents a minimal feed- 
back model with its main properties and results. Prelim- 
inary results were presented by Monaco (2002; 2003). An 
upcoming paper will focus on the destruction of collapsing, 
star- forming clouds (Monaco 2004, hereafter paper II). 

The paper is organized as follows. Section 2 describes 
the physical ingredients of the model. Section 3 introduces 
the system of equations used. Section 4 the main solutions. 
Section 5 is devoted to a discussion of the results, and Sec- 
tion 6 gives the conclusions. Finally, three appendices give a 
list of frequently used symbols, a determination of the time 
scales of coagulation of cold clouds and a study of the fate 
of SBs in the Uh-ias plane. 



2 FEEDBACK BY STEPS 

Feedback is assumed to take place through a chain of pro- 
cesses: 

(i) The densities and filling factors of the two phases 
are determined by pressure equilibrium. 

(ii) The cooled or infalled gas fragments into clouds 
with a given mass spectrum; this is truncated at low masses 
(which are easily destroyed) and at high masses (which con- 
tinually collapse). 

(iii) Collapse is triggered in clouds larger than the Jeans 
mass; we use a criterion valid for non-spherical clouds. 

(iv) Collapsing clouds are continually created by coag- 
ulation. 

(v) Stars form in collapsing clouds. Self-regulation of 
star formation by HII regions destroys the clouds before 
most SNe explode. 

(vi) SN remnants (hereafter SNRs) soon percolate into 
a SB, which swoops the ISM. SBs heat the gas whenever they 
are in the adiabatic stage, i.e. when the interior gas has not 
had time to cool, while they collapse (and thus cool) the 
hot phase into a thin cold shell whenever they get into the 
so-called snowplow stage. 

(vii) SBs stop sweeping or collapsing the hot phase when 
they remain pressure-confined or overtake the typical verti- 
cal scale- height of the system (blow-out). 
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In the following we describe these steps in detail. All 
distances are given in pc, masses in M©, times in yr, tem- 
peratures in K, gas densities in cm~^, average densities in 
Mq pc~^, surface densities in Mq pc~^, energies in lO^"' erg, 
mechanical luminosities in 10'^** erg s^^, mass flows in M0 
yr"'^, energy flows in 10^^ erg/yr. Pressures are divided by 
the Boltzmann constant k and given in K cm~^. 



2.1 Pressure equilibrium 

Let's consider a volume V filled with a two-phase medium, 
with temperatures of hot and cold phases Th and Tc and 
densities rih and ric- The volume is assumed to be large 
enough to contain many star-forming clouds. An external 
halo acts as a reservoir of gas, which continually replenishes 
the cold component^. Stars form from the cold gas. The four 
components (cold and hot phases, stars and the external 
halo) have masses Mcoid, Mhot, and Mhaio- The total 
mass of the system is fixed to Mtot- The temperature of the 
cold gas is kept fixed to 100 K, i.e. roughly the position 
where the cooling function of the gas drops, so that further 
cooling is inhibited unless the cloud collapses and its density 
gets very high. Let jUh and /Xc be the mean molecular weights 
of the two phases, /h and /c their filling factors (/h + /c = 1), 
Ph = Mhot/V and pc = M^oid/V their average densities and 
Fh = A/hot /(Afcoid + Afhot) the fraction of hot gas. Pressure 
equilibrium implies: 



From this we obtain: 
f = I 



(1) 



(2) 



and of course /h = 1 — /c, = ph//hp.h"i-p and ric = 
nc//cA*cWp (where mp is the proton mass). Finally, the de- 
pendence of the Ph and pc molecular weights on metallicity 
is taken into account. 



2.2 Fragmentation of the cold phase 

It is assumed that the cold phase fragments into clouds with 
a given mass spectrum. As commented in the introduction, 
according to the turbulent picture of the ISM the "clouds" 
(i.e. peaks of the fractal density fields) are not stable entities 
but transient features of the medium. We will assume in the 
following that the self-gravitating clouds are reasonably sta- 
ble (in the sense that they are not significantly reshuffled by 
turbulence) within one or two dynamical times and that the 
continuous reshuffling of the density field does not change 
the statistics of clouds. 

The mass spectrum of the so-defined clouds is assumed 
to be a power-law: 



Na{mci)dm^i = Afo(mci/l M©) ""'dmd : 



(3) 



where A^o is a normalization constant (with dimensions pc ^ 
Mq^), fixed by requiring = J Ncimcidrrici (see below), 



^ The halo is assumed, for simplicity, to be completely decou- 
pled from the hot phase, although in reaUstic situations the two 
components will interact. 



and aci is a free parameter. This choice is the natural out- 
come of many different processes, including turbulence. The 
parameter ad can be constrained both from theory and ob- 
servations of the ISM (see, e.g., Solomon et al. 1987), and 
should vary between 1.5 and 2 (the latter considered as a 
reference value), at least in self-regulated situations like the 
Milky Way. Notice that in this way a significant amount of 
mass is located in high-mass clouds. 

To the clouds we associate a typical radius ad defined 
simply as ma = 47rOdPc/3, or: 



md = 0.104 pcTicad M© . 



(4) 



This does not imply an assumption of sphericity of the 
clouds. 

The mass function of clouds is truncated both at low 
and high mgisses. At the high mass end the mass func- 
tion is truncated by gravitational collapse, because clouds 
that form stars are quickly destroyed. The upper mass limit 
rUu will be computed in the next session. At low masses 
clouds are easily destroyed by a number of possible pro- 
cesses, among which thermo- and photo-evaporation. Mc- 
Kee & Ostriker (1977) set the lower limit to Od = ai = 0.5 
pc. For pc~1.2 and nc~10 cm~^ this corresponds to mi~0.1 
Mq . We set the lower mass limit to this value. This is surely 
a rough approximation, as mi should be self-consistently de- 
termined by the dynamics of the system, and is unlikely to be 
a constant. However, its actual value does not have a strong 
impact on the results as long as rriu mi, a condition that 
is verified by most solutions. Nonetheless it is important to 
set mi to a non-vanishing value both to avoid divergence in 
a few calculations (like the normalization of the mass func- 
tion for ad^2) and to avoid contributions from clouds that 
most likely do not exist. 

The normalization constant of the mass function is: 

Pc 



Ar„(l Mo)"-^' = 



/(mu,mi) 



(5) 



Here the function /(mu,mi) is equal to (mu""'"*"^ 



)/(— Od + 2) if Qci 7^ 2, otherwise /(mu,mi) 



ln(mu/mi). 



2.3 Critical mass for clouds 

Massive clouds are destroyed by gravitational collapse. In 
absence of magnetic fields and turbulence the threshold mass 
for collapse is fixed by the Bonnor-Ebert criterion (Bonnor 
1956; Ebert 1955), and depends on an external pressure term 
Pcxt- If the external pressure is fixed to the thermal one, 
the criterion is equivalent to the classical Jeans mass. To 
generalize it to non-spherical clouds, we follow Lombardi & 
Bertin (2001), who find: 



mj 



1.18- 



(6) 



Here Cs,c is the sound speed of the cold phase, the external 
pressure is set to the thermal one and the parameter pshape 
is defined by the authors as: 



Pshape ■ 



127r 



3 N 1/3 V^4/3 



52 



(477) 



(7) 
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( I |V5M(sx)r'dS / |Vjw(sx)|ds] 
\Jav Jdv J 



In this equation the integrals are performed on the sur- 
face dV (of area S) of the volume V of the cloud; the 
function u is the cloud density normalized to its maximum 
value u = p/pmax and is a "Jeans length" defined as 
s = 1^ 47rGpmax/ ci,c ; X is the space coordinate and ^ = sx. 
The parameter /ishapo is dimensionless, scale invariant (i.e. 
does not change for similarity transformations) and is always 
smaller than unity. For a sphere jUshape = 1, and the Jeans 
(Bonnor-Ebert) criterion is recovered. In general, collapsing 
clouds will be non-spherical, and this will correspond to an 
increase of the threshold mass mj. We treat /ishapc as a free 
parameter. It can be considered as a product of two terms, 
/ii and ij,2 , given in the first and second lines of Equation 7. 
Both terms are < 1 and are unity for a sphere; moreover, fj,2 
is unity when gravity is negligible. So, a rough estimate can 
be obtained by considering /ishape~pi. 

This quantity is easily computed in the simple case of a 
rotational ellipsoid with semi-axes ai and a2 (with the third 
semi-axis 03 = 02). If r = a2/oi is the axial ratio, we find: 



Mshape 



t?(r)2r4/3 ' 



(8) 



where g{r) — 1/2 + arcsin \/l — r'^/2r\/l — if r < 1 and 
g{r) = 1/2 + log[(r + Vr^ - l)/(r - Vr^ - l)]/4rVr2 - 1 if 
r > 1. In this case /Xshape takes values ~0.5 for axial ratios 
of order 1:5 (in both senses), while it gets to ~0.2 for axial 
ratios 1:10. As this is likely to be an overestimate of the 
actual value, we consider 0.2 as a reference value for this 
parameter. 

Magnetic fields and turbulence could in principle inval- 
idate the Bonnor-Ebert criterion by providing non-thermal 
support to the cloud. Recent simulations (see, e.g., Mac Low 
2003) have shown that turbulence cannot inhibit the collapse 
of critical clouds; the Jeans criterion remains valid provided 
that the quadratic sum of kinetic and sound speeds is used in 
place of the sound speed itself. For a typical turbulent speed 
of several km s^^, the Jeans mass would correspond to that 
relative to a temperature Tc of several 10^ K. The effect of 
turbulent motions can thus be roughly implemented by as- 
suming a very small value for yUghape, of order 0.01. Magnetic 
fields can halt the global collapse of the cloud but not its 
fragmentation into stars, so their effect on the critical mass 
for collapse is negligible. 

Finally, in cases like the sweeping of a spiral arm or 
during a merger the Jeans criterion can be changed by ex- 
plicitly introducing a Pcxt term. This will correspond to a 
sudden decrease of the Jeans mass, and then to a burst of 
star formation. 



2.4 Coagulation of cold clouds 

Clouds larger than the Jeans mass are continually created 
by kinetic aggregation (coagulation) of smaller clouds. This 
is described with the aid of the Smoluchowski equation (von 
Smolucliowski 1916). In this we follow the approach of Cav- 
aliere, Colafrancesco & Menci (1991; 1992; see also Menci 
et al. 2002), who used this formalism to describe the kinetic 
aggregation of dark-matter halos. 



The details of the calculations are reported in Appendix 
B. In brief, the coagulation of clouds is driven by a kernel: 

K = pc ((EcoagWap)^)^ . (9) 

Here Ecoag is the cross-section for interaction and Uap is the 
approach velocity, while the two averages are done over ve- 
locity and mass. Notably, it is assumed that clouds, although 
transient, are stable for one crossing time ad^ap; this is rea- 
sonable as Vap is typically larger than the sound speed of the 
cold phase. Following Saslaw (1985) the cross-section for the 
coagulation of two clouds (denoted by 1 and 2) is: 



n{ai+a2)^ 1 -|- 2G 



, (mi + m2) 



tti 



02 



(10) 



The first term corresponds to geometric interactions, the 
second to resonant ones; this last term is effective when the 
approach velocity is not much larger than the internal veloc- 
ity dispersion of the clouds. In most cases considered here 
the geometrical term results dominant, so we will neglect 
resonant interactions in the following. Notice that this cross- 
section is valid for spherical clouds; we do not consider the 
effect of asphericity here, as it would be a further-order cor- 
rection with respect to that of the Jeans mass introduced 
above. 

It is shown in Appendix B that the time scale for coag- 
ulation is: 



Pc (Uap) 



_ (Any'-' 1 

tcoag - [-jj -P. 



(11) 



The typical mass scale of the mass function, identified with 
the upper cutoff, grows like (l-|-i/3tcoag)^. For a Maxwellian 
distribution of velocities with ID dispersion Uv we have 

(Vap) = l.SOcTv. 

The time at disposal for accretion is the time necessary 
to a Jeans-mass cloud to be destroyed. This will be related 

to the dynamical time: 



tdy 



32Gpc 



;5.15x 10^ (/Ucnc)"^''^ yr. 



(12) 



As star formation is triggered roughly after tdyn, and eaxly 
feedback from young stars destroys the cloud in a compara- 
ble time (see below), we conservatively allow aggregation to 
go on for two dynamical times. Thus, the upper mass cutoff 
is set to: 



rriu = mj 1 -h 



2t. 



dyn 



The mass of the typical collapsing cloud is then: 
XTj" i^ciNc\{md)dmc\ 
XT" Nci{mci)dmci 



(13) 



(14) 



and the fraction of cold gas presently available for star for- 
mation is: 



fa 



f"" mciArci(mci)dmci 



The total number of collapsing clouds is: 

-/Vcc = /coll 



(15) 



(16) 
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Coagulation is a physically motivated and reasonable 
mechanism to explain the growth of cold clouds, but it has 
never been validated (to the best of our knowledge) by sim- 
ulations that include MHD turbulence. Besides, it heis been 
proposed that giant molecular clouds form in the converging 
flows caused by the sweeping of spiral arms (Ballesteros- 
Paredes, Vaaquez-Semadeni & Scalo 1999), a process that 
cannot be introduced without a proper modeling of the disc. 

A consequence of the assumptions done is that cool- 
ing alone is not going to produce clouds larger than the 
Jeans mass; they are produced only by coagulation of smaller 
clouds. This is contrary to the naive expectancy of a mass 
function of clouds which is truncated below by the Jeans 
mass in case of a cooling flow, as only fluctuations larger 
than the Jeans mass can grow. This is not what is observed 
in the case of thermal instability in turbulent media, where 
(without thermal conduction and UV heating) structures of 
all masses arc observed down to the resolution limit (see, o.g, 
Kritsuk & Norman 2002). On the other hand, it is possible 
that giant clouds, much larger than the Jeans mass, form 
in the cooling flows that take place at the centres of cosmo- 
logical halos. This is neglected here, but can be modeled by 
introducing a further mass scale in the mass function. 

2.5 Star formation and early feedback 

Collapsing clouds can reach high enough densities to trig- 
ger the formation of H2 and further cool to ~10 K. After 
one dynamical time (Equation 12) star formation starts in- 
side the "molecular" cloud. An important point is that early 
feedback from massive stars can destroy the collapsed cloud 
before the bulk of type II SNc has exploded. It has been 
shown (Franco, Shore & Tcnorio-Tagle 1994; Williams & 
McKee 1997; Matzner 2002) that HII regions are a source of 
turbulence, and their energy input is sufficient to destroy the 
star-forming clouds, pre-heating them at >10'' K. A similar 
role is played by stellar winds, that are typically trapped 
inside HII regions (McKee, van Buren & Lazareff 1984). 
Matzner (2000) computed the amount of turbulence driven 
into the star-forming cloud by expanding HII regions. Un- 
der the assumption that the rate of injection of turbulence 
equates the decay rate estimated from N-body simulations, 
he predicted that the cloud would be destroyed in ~2 x 10^ 
yr, i.e. about one dynamical time of the uncoUapsed cloud 
(Equation 12 with nc~10 cm~^), with a resulting efficiency 
of star formation /, (i.e. the fraction of the cloud that goes 
into formed stars) of ~5-10 per cent. This is in rough agree- 
ment with both observations of molecular clouds and esti- 
mates from globular clusters (/*~1-10 per cent; see, e.g., 
Elmegreen 2000, 2002). 

The ability of the energy from SNe to emerge from the 
destroyed cloud, possibly the most delicate step in the whole 
chain of feedback events, is addressed in paper II; here we 
give only a very short summary of the results. When SNc 
start to explode the cloud is already in the process of being 
destroyed, so that a significant fraction of mass is in a warm, 
diflfuse phase. SNRs propagating in this dense environment 
soon radiate their thermal energy (see next section for more 
details). In this case, the mass internal to the blast collapses 
into a thin, dense shell that fragments as soon as the blast 
is confined by kinetic pressure. So, the net effect of the first 
SNe is that of collapsing again the diffuse material heated 



up by the HII regions. After a few SNc, most gas is re- 
collapsed into cold clouds with a low filling factor, while the 
diffuse component has such a low density that SNRs emerge 
from the cloud before cooling. From this point all the energy 
from SNe is used to drive the SB. In case many tens of SNe 
explode in a single cloud, most energy (90-95 per cent) is 
used to drive the SB, while for very small clouds, where only 
a few SNe explode, the first SN is able to destroy the cloud, 
losing most of its energy in the process, while the other SNe 
(if any) will pump energy into the ISM with a likely high 
efficiency. Eventually, only ~10 per cent of the initial cloud 
is found in diffuse, hot gas with temperature of order 10^ K; 
lower values arc expected if the cloud is particularly dense. 

For this version of the feedback model we decide to give 
a minimal, heuristic description of this process, in order to 
keep the model simpler. Each SN releases lO^^Ssi erg in 
the ISM^ We assume that all the energy is available for 
driving the SB; in case of very small collapsing clouds a 
lower effective value of Er,i will be plausible. We assume that 
a fraction /evap of the cloud is evaporated to a temperature 
Teva.p, while the rest (amounting to a fraction 1 — /* — /evap) 
is re-collapsed into cold clouds. Of course /evap+/*<l. We 
use as reference values £51 = 1, /cvap=0.1 and rcvap=10^ K, 
with the warning that in case of very dense clouds /evap will 
likely be lower (see paper II). 

Finally, the contribution of a single collapsing cloud to 
the global star formation rate is: 



mcc 

^dyn 



(17) 



2.6 Super-bubbles 

SNRs associated to massive stars in a star-forming cloud will 
soon percolate into a single hot bubble. As a consequence, 
all the SNe exploding in a cloud will drive a single SB into 
the ISM (see, e.g., Mac Low & McCray 1988). 

Stars are formed with a given Initial Mass Function 
(hereafter IMF) that must be specified. For the model the 
only information needed is the mass of stars formed for each 
supernova, M*,sn. We associate one SN to each > 8 Mq 
star; if the (differential) IMF has a slope — (aimf + 1) and 
the lifetime of a star goes like its mass raised to — aufc, the 
rate of SN explosion goes like f("imf-«iite)/aiito p^j. standard 
choices of aimf = 1.35 and aufe ~ 2.5 — 3 the exponent takes 
a value of ~ —0.5. In other words, the rate of SN explosion 
depends weakly on time, and is approximated as constant. 
Denoting by tufa the difference between the lifetime of an 
8 Mq star and that of the largest star, the number of SNe 
that explode in a collapsing cloud and the resulting rate are: 



TVs, 



/* 



(18) 



(19) 



tlife-^*,sn 

The mechanical luminosity of the SB is then Lmech = iss x 
10"^* erg s^^, where: 



^38 = 



10"j?sng51 

1 yr 



(20) 



Observations suggest values of E51 in the range 1 to 10. 
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Adiabatic stage 



-Tit) 

Post-shock temp. T^/^'*' (i) 



Radius 
Shock speed 
Average temp. 



Post-shock press. 
Coohng time 
Swept mass 



cool ^ ' 



81.3 {L38/f^i,ni,y/Hl^'^ pc 
47.7 {L38^h«h)^/'^*6 km s-i 
1.79 X 10'^' L%'l^'J'n-'/%'/' K 

-2/5 ,—• 



5.4 X 10* V^^^nr^/^ij'/' K 



2.16 X 10^ Ll(^^{in,ni,y^/->t^ ^'^ K < 

255 7-3/5 9/10 -8/5,-1 -6/5 

5.53 X lO'' L]l^{ixwn^fl^r'J^ M© 



Internal mass 


Mint(t) 


= Msw 


PDS stage 


Radius 




= 70.2 (L38/Mhnh)i/5t6^^ pc 


Shock speed 




= 41.2 iL;i8/fii,ni,)^/H-^^^ km s"! 


Post-shock press. 




= 1.60 X 10^ Ll{^^{fn,ni,f/^t~'^^^ K cm-3 


Swept mass 


Msw(t) 


= 3.56 X lO* L|/'^(/ihnh)2/5t^/^ Mg 


Internal mass 




= Msw (1 - (Ws)"^-^) 



Table 1. Main properties of SBs. 



In presence of a two-phase medium the SB expands into 
the more diffuse, more pervasive hot phase; cold clouds will 
pierce the blast, but this will promptly reform after the cloud 
has been overtaken (McKee & Ostriker 1977; Ostriker & 
McKee 1988; Mac Low & McCray 1988). 

The evolution of the SB is described following the model 
of Weaver et al. (1977; see also Ostriker & McKee 1988). In 
the beginning the SB is adiabatic, because the shocked ISM 
has not had time to cool. In this case: 

i^r(^)=8l.3fi^y'^^pc, (21) 

where tfi = yr. Table 1 reports the main properties of 

the SB expanding in the hot phase. 

In the adiabatic stage the hot phase is shock-heated by 
the blast. Of the initial energy of the SN, 73.7 per cent is 
thermal and 26.3 per cent is kinetic. This stage ends when 
the post-shock mass elements cool, i.e. when t^ooi(*) ~ 
The cooling time is computed as: 

tcooi = 3fcr/nhA(T) (22) 

and is evaluated at T^^'^^ (given in Table 1) and 4nh (due 
to the shock jump condition). For the cooling function we 
use the approximation proposed, e.g., by CiofR, McKee & 

Bertschingcr (1988): 

A = 1.6 X lO-^'CmT"'/', (23) 

where ("m = Zhot/ZQ is the metallicity of the hot gas in 
solar units. This formula is relatively accurate in the range 
more realistic cooling function would 
be desirable, but would make analytic estimates unfeasible. 
The time of shell formation is then: 

tpd. = 2.33 X 10*L|/"M^/^^n-«/"C-^/" yr (24) 

We call Rpds the radius of the SB at tpda- After this moment 
the swept mass collapses into a thin cold shell. This shell acts 
like a snowplow, making the swept ISM collapse into it. For 
simplicity the cold clouds are assumed as before to pierce 
the shell without any effect. Some of the hot gas will any- 
way remain inside the bubble, pushing the snowplow with 



its pressure; this stage is called Pressure Driven Snowplow 
(PDS). We use the solution of Weaver et al. (1977; see also 

Castor, McCray & Weaver 1975) that includes thermal con- 
duction at the interface between the gas and the cold shell, 
a mechanisms that releases more hot gas from the shell into 
the interior. They obtain: 

Ril^'Ht) = 70.2 ( ' tl^'pc, (25) 

Notice that the time dependence is the same as above, due 
to the presence of an increasing amount of hot interior gas. 
This gas is however negligible with respect to the swept 
mass, and is so diluted that further cooling is inhibited. A 
more standard choice for the evolution of the SB in the PDS 
stage would be (see, e.g., Koo & McKee 1992) 7?^^'*"^ oc t'^^''; 
the exponent decreases only by 5 per cent with respect to 
Weaver et al. (1977). Table 1 reports the main characteris- 
tics of the SBs in this stage. 

To ease numerical integration we interpolate between 
the adiabatic and PDS stages assuming that after tpds the 
blast radius evolves like Rsh oc t"'^ and the velocity like 
Vsb OC until the PDS solutions are met. 

In the PDS stage, the amount of ISM swept by the SB 
that is collapsed into the shell is estimated as the fraction 
of the internal material for which (in the adiabatic solution) 
icooi('";^) < Assuming a power-law profile for density and 
temperature of the gas just inside the adiabatic blast, in 
the pressure-gradient approximation of Ostriker & McKee 
(1988) we obtain T oc {r/R^b)"-^ and p oc {r/R^b f. From 
these relations we obtain that the internal mass (not yet 
collapsed into the shell) is related to the swept mass as: 

Mint = Msw (^/^pdB)"^■^ (26) 

This is valid of course for t > ipds. For simplicity we assume 
that the thermal energy of the SB is lost at the same rate: 

E^^^ = 0.737 Rsnt (l - (t/tpd,)-''-^) , (27) 

while the kinetic energy is kept at E^l''''^ = 0.263i?snt. 

The explosion of the last SN marks the exhaustion of 
energy injection into the SB, so the evolution after this 



Physical regimes for feedback 7 



event should follow that of a SNR. We observe that SBs 
are stopped by thermal pressure or by blow-out (see below) 
before exhaustion in virtually all cases, so an accurate mod- 
eling of this stage is immaterial. In any case, we assume that 
after the last SN has exploded the blast always evolves like 
the adiabatic Sedov solution for a SNR, Rsh oc t^^^ . 

A note of caution is necessary on the application of 
these solutions for the evolution of the SB. They are valid 
if the hot phase is uniform, the cold phase negligible, and 
if the mass of the "wind" that drives the SB is negligible 
with respect to the swept mass. This last condition is vio- 
lated in most actual cases as soon as a significant fraction of 
the collapsing cloud is evaporated. On the other hand, the 
other conditions arc also violated: the ISM is structured, and 
events like thermo-evaporation of clouds, dragging of clouds 
by the internal gas, turbulent and magnetic pressure and 
cosmic rays are likely to influence significantly the dynam- 
ics of the SB. Some of those effects can be modeled analyti- 
cally (see Ostriker & McKcc 1988), but at the cost of more 
uncertainties and dangerous assumptions. We decide to rely 
on the simple solutions given above, with the caveat that all 
numbers must be considered as useful order-of-magnitude 
estimates. 

Nonetheless, we have tried to include thermo- 
evaporation of the cold phase by the expanding SB, so as to 
quantify the mass flow implied. Following the approach of 
McKee & Ostriker (1977) and the generalization of Ostriker 
& McKee (1988) to SBs, wo find that the thcrmo-cvaporatcd 
ISM is generally much less massive than the evaporated gas 
from the star-forming cloud, so the inclusion of thermo- 
evaporation, while introducing further uncertainties, does 
not affect strongly the results presented here. We will ne- 
glect thermo-evaporation in the following. 

2.7 The fate of SBs 

An SB can end in two ways: (i) being confined by external 
pressure, (ii) blowing out of the system. 

Case (i) takes place at time tconf when the shock speed 
is equal to the external, thermal one: 

Wsb(*conf) = Cs,h = 91.2(Th,6/Ath)'''^ (28) 

where Th,6 = Th/lO^ K. As the blast propagates into the 
low-density hot phase, kinetic pressure is always negligible. 
The time and radius at which confinement takes place are 
given in Table CI of Appendix C. After confinement, the 
blast (in the adiabatic stage) dissolves in the hot phase or 
the shell (in the PDS stage) fragments because of Raileigh- 
Taylor instabilities. This allows the hot phase to mix with 
the interior hot gas. However, as long as tconf < tiifo many 
SNe explode after confinement. This will correspond to the 
creation of secondary bubbles; the medium in which they 
expand will depend on the velocity with which the inte- 
rior gas mixes with the external one. In the adiabatic con- 
finement case, it is easy to see that secondary bubbles will 
be confined in the adiabatic stage as well, so all the en- 
ergy from SNe will be released to the hot phase; in this 
case feedback is mostly efficient. In case of confinement in 
the PDS stage, the situation is more complicated. If inte- 
rior and external hot gas mix very quickly, the secondary 
bubbles will expand in the same medium and will then cre- 
ate secondary shells, but if mixing is slow the energy of the 



remaining SNe will be pumped efficiently into the hot, rar- 
efied internal gas. To address this case we assume that the 
energy of the SNe exploding after confinement is released 
to the hot phase either entirely, /pds=l, or by a fraction 
/pds=0. 737(1 - (tconf /tpds)"'''-^) + 0.263, that takes into ac- 
count that thermal energy is dissipated according to Equa- 
tion 27. The two cases should bracket the true solution. 

Case (ii), the blow-out of the SB, takes place when the 
SB overtakes the vertical scale-height Hf,s of the system, 
defined as (Mac Low & McCray 1988; Koo & McKee 1992): 

H^ff= — / pi,{z)dz, (29) 
Po Jo 

where z is the vertical direction (that for which Hgfc is min- 
imal) and po = ph{z = 0). The blow-out condition is obvi- 
ously -Rsb = HcB, and the blow-out time is reported in Table 
CI. This condition is true if all SBs are centred at z = 0. 
Coagulation naturally leads to mass segregation, and this is 
in line with the observational evidence that molecular clouds 
show a smaller vertical scale-length then HI. However, typ- 
ical blowing-out SBs will be away from the mid-plane and 
will blow out only from one side, with a result that is not 
vastly different from a bi-polar blow-out. Moreover, off-plane 
SBs will blow out more easily, and intermediate configura- 
tions with mid-plane SBs being pressure-confined and exter- 
nal SBs blowing out will be possible. In order to keep the 
model as simple as possible we consider only mid-plane SBs; 
the Hes vertical scale-length will then be understood as the 
difference between the HI and H2 scale-lengths 

The SB does not stop immediately after blow-out, as 
the rarefaction wave that follows blow-out must have time 
to reach the blast traveling in the horizontal direction. We 
then allow for a sound crossing time before stopping the 
SB. If is the mean effective Mach number of the blast 
(the square ratio between the blast speed and the average 
internal sound speed) and Rsh{t) oc f, the sound crossing 
time of an adiabatic bubble is: 
cx 

Across — ^bo j (^0) 

V 

For adiabatic and PDS blasts, = 1.61 and 1.18 (Ostriker 
& McKee 1988; Weaver et al. 1977). The final radius R^o' 
and time tbo' (= ^bo + Across) at blow-out are reported in 
table CI. Between tho and tbo' the SB can get into the PDS 
stage (if it hasn't yet) or be confined by pressure. 

The final time tfm and radius Ran are defined respec- 
tively as the smallest between ty,o' and tconf, and between 
J?bo' and /Jconf. Appendix C reports a study of the final 
state of SBs (confinement or blow-out in adiabatic or PDS 
stage) in the plane defined by the two variables Uh and Lag. 

At blow-out part of the hot interior gas of the SB es- 
capes to the halo. To compute the fraction of hot blown- 
out gas we adopt the following simple geometrical model 
(Fig. 1). The ISM swept by a SB of radius _Rfin blows out 
from the two polar cups defined by the intersection of the 
sphere and the two horizontal planes at distance HeS from 
the centre. The swept gas receives momentum from the blast 
in the radial direction, so the blowing-out gas is that con- 
tained in a double cone with the opening angle 6 of the 
polar cups; we have cos 6 = Hes/Rfin- Neglecting the ISM 
contained in the polar cups (that are outside the volume V), 
the fraction of swept ISM that is blown out is: 
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Figure 1. Geometrical model for blow-out. The SB starts blow- 
ing out when its radius is equal to -ffcffi but continues to expand 
for one sound crossing time, finally reaching the radius -Rfln. The 
two polar cups (diagonal-shaded regions), defined by the inter- 
section of the final SB and the two planes at -ffcffi are assumed 
to be devoid of matter. All the matter present in the double cone 
(vertical-shaded regions) with an aperture equal to that of the 
polar cups receives a radial momentum that allows it to blow-out 
into the halo. 

f _/ \[H^s/Riir.-{H,a/Rfir.f] if Rsn > H^a 

^"""'X if i?fln<i/cff. ' 

This is valid both for adiabatic and PDS blow-out. Consis- 
tently, the absence of ISM in the two polar cups is considered 
when computing the swept mass. With this simple model, 
that contains no free parameters, the fraction /bo ranges 
from to ~0.2; this is roughly consistent with Mac Low & 
McCray (1988) and Mac Low, McCray & Norman (1989), 
who report that most of the internal hot gas remains in the 
disc. It is anyway interesting to allow for higher /bo values; 
this is done by forcing the maximum of Equation 31 to be 
/bo, max, which is takeu as a free parameter. 

A note on definitions; in a cosmological context the term 
blow-out is used for the gas that is expelled from a galactic 
halo; we use it for the expulsion of gas from the "galaxy", 
i.e. from the region where stars and ISM are present, but our 
blown-out gas is destined by construction to remain in the 
halo. As already mentioned above, this over-simplification 
is introduced to avoid modeling of the external halo. As we 
know temperature, density and escape velocity of the blown- 
out gas, modeling of galactic winds is readily feasible once 
the global properties of the hosting dark-matter halo are 
specified. 

The energy of the SNe exploding after tbo' is assumed 
to be funneled out into the halo. Nonetheless, the restored 
mass, responsible for chemical enrichment, is for simplic- 
ity blown out with the same efficiency (/bo) as the rest of 
the mass, leading to a possible underestimate of the metals 
ejected into the halo. We will see in the following that the 




Figure 2. Mass flows between the four components described 
in the model. Arrows denote the flows connected to infall (Mj^f), 
star formation (Mgf), restoration (Mrcst), cooling (McooOi evapo- 
ration (Movap), snowplows (Mgnpi), leak-out (Mioak) a^nd the rate 
at which the hot phase is engulfed by SBs (Mint). Blow-out takes 
mass by a fraction /bo from the internal, evaporation, snowplow 
and restoration mass flows. 

ejection of metals in the halo is anyway rather efficient even 
with this assumption (see also de Young & Gallagher 1990; 
Ferrara & Tolstoj 2000). We leave a refinement on this point 
to further work. 

The porosity Qsb of the SBs is defined as the fraction of 
volume occupied by expanding blasts. This is a very impor- 
tant quantity, as when it is unity it indicates that the blasts 
percolate the volume and create a super-SB. The computa- 
tion of Qsh depends on the time at which SBs stop to expand 
into the ISM. In case of blow-out the blast halts at tbo' and 
the energy of further SNe is funneled out of the volume V 
into the halo. On the contrary, in case of confinement sec- 
ondary blasts form after tan, whose energy is still injected 
into the ISM. To take this into account we compute Qsb as: 

0.b = -^-7r/ RL{t)dt. (32) 

tdynV 3 /„ 

where tpao ~ tfin in case of blow-out, or fufc in case of 
confinement; in the latter case Rah is kept constant to Ran 
after tfi^- Regarding the adiabatic blow-out regime, it is clear 
that if SBs remain identifiable after fporo, then a value of Qsh 
inferred from observations will be higher than that given by 
Equation 32. However, recognizable bubbles do not play the 
same dynamical role as expanding blasts. 



3 THE SYSTEM OF EQUATIONS 
3.1 Mass flows 

Fig. 2 shows all the mass flows between the four components 
that are taken into account in this model. 
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Cold gcis is continually infalling from the halo. This is 
modeled simply as follows: 

Mhalo 



Mint = 



tint 



(33) 



where tin! is a parameter of the system. 

The hot phase cools at the rate tcooi; if it were com- 
pletely homogeneous it would remain globally hot. In real- 
istic cases the hot phase will show a rather broad range of 
densities and local cooling times, so that a fraction of the 
gas will be able to cool to low temperature. A modeling of 
this fraction would require detailed knowledge of the density 
distribution of the hot phcise; we prefer to leave it as a free 
parameter, /cool, with 0.1 as a tentative reference value. The 
cooling mass flux is then: 

Mhot 



= /c< 



(34) 



^cool 

where tcooi(7h) is computed from Equation 22 and with the 
simple cooling function given by Equation 23. 

While the cold phase is easily confined by a modest 
gravitational well, the hot phase is generally able to leak 
out of the volume V to the halo. The time scale connected 
to this leak-out is: 



tleak = Y 3/d 



Cs,h 



(35) 



where Ca,h is the sound speed of the hot phase and d is unity 
if leak-out is on one preferential direction, 3 if it is spherically 
symmetric. In the following we will consider for simplicity 
leak-out in one direction; we have verified that the results 
are not sensitive to the value of d. The mass loss rate is then: 

Mhot 



Ml 



^leak 



(36) 



This term should be revised if some external hot halo gas 
hampers leak-out. 

Cool gas transforms into stars at the rate: 



Msf = /*/coll 



M, 



cold 



tdy 



(37) 



This is easily obtained by multiplying the contribution of a 
single cloud (Equation 17) by the total number of collapsing 
clouds A'cc (Equation 16). It can be written also as Mst = 
Mco\d/ta!, where 



^sf — idyn/ /*/coll - 



(38) 



A fraction /rest is instantaneously restored to the hot 
phase: 



Mr, 



f rest Mat . 



(39) 



This flux is responsible for chemical enrichment; we notice 
that this equation implies instantaneous recycling. 

The rate at which the mass of collapsing clouds is evap- 
orated back to the hot phase is: 

Afcold 



M, 



fcvapfcoU ' 



dyn 



(40) 



It follows that Mevap = feva.pMsf/fi,. 

At the final time each SB has swept a mass 
Msw(tfin), of which a part Afint(tHn) (sec Table 1 and Equa- 
tion 26) is in hot internal gas and the rest is in the snowplow. 
The rate at which the hot phase becomes internal mass of a 
SB is: 



iVc 



Mint (fan) 
^dyn 



while the rate at which it gets into a snowplow is: 

Msw(fan) - Mint (ten) 



Msnpl = Nc, 



(41) 



(42) 



We recall that a fraction /bo (Equation 31) of the 
swept material (both hot and cooled) and of the restored 
and evaporated mass is blown-out to the halo. Defining 
Mbo = /bo (Mevap + Mrest + Mnt "I" Msnpi), the System of 
equations that describes the mass flows is: 



( Mcold = Mint + Mcool - Maf - Mevap 
+ {1 - /bo)Afsnpl 
Mhot = -Mcool - Msnpl - Mleak - /boMnt 

+ (1 - /bo)(Mevap + Mrest) 
M* = Msf- Mrest 

I Mhalo = -Minf -h Mleak + Mbo 



(43) 



Mass conservation is guaranteed by the condition Mhot + 

Mcold + M* + Mhalo = 0. 



3.2 Energy flows 

A similar set of equations can be written for the energy 
flows. Here we concentrate only on the energy Shot of the hot 
component, that determines Th. The total energy releEised 
by SNe is: 



tdyn 



(44) 



The rates of energy loss by cooling, snowplow, blow-out and 
leak-out are respectively: 



■E'snpl = MsnplTh 



-Ebo — /boMintTh 



k 



2 Hhmp ' 

3 k 
2 Hhmp ' 



Sleak 



-Ehot 
^leak 



(45) 



(46) 



(47) 



(48) 



Regarding the energy budget of the SB (Equation 27), 
while thermal energy is obviously given to the ISM kinetic 
energy is transformed into turbulence and then partially 
thermalized. We have verified that including or excluding 
kinetic energy from the energy budget does not change ap- 
preciably the dynamics of the system. At present we decide 
to give it to the ISM, more refined modeling of the decay of 
turbulence will be presented elsewhere. 

In the blow-out regime the ISM receives a fraction (1 — 
/bo) of the energy of the SB (thermal. Equation 27 evaluated 
at fan^, plus kinetic) and of the energy connected to the 
evaporated mass (that comes from the first SNe exploding 
in the cloud): 



Efb = (l-/bo) K 



E. 



(th) 



I 77, kl 



(kin) 



tdy 



\ -^i^-'evap-'-evap c\ 

2 AAhmp 



.(49) 



When ten is short we force iJsnten to be at least unity. 
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In the adiabatic confinement case all energy from SNe is 
given to the ISM: 

Sfb = 4n . (50) 

In this case the energy of the evaporated cloud is already 
included in the total SN budget. Finally, in case of PDS 
confinement the ISM receives the thermal and kinetic energy 

of the SB. As mentioned in Section 2.7, the energy -Brest 
associated to SNe exploding after is given either by a 
fraction /pdB=0.737(l - (tflnApds)"^ '') + 0-263 (the case of 
fast mixing, where secondary bubbles share the same fate as 
the principal one; see Section 2.6) or entirely, /pds=l (the 
case of slow mixing, where the remaining energy is pumped 
into the rarefied interior of the bubble): 

Bfb = (sir' + -^ib"' + /pdsSrest) (51) 

t'dyn 

Given the uncertainty connected to the modeling of -Brest, 
we consider it worthless to include a detailed treatment of 
the energy from the evaporated cloud in this case. 
The equation for the evolution of Bhot is: 

B'hot = — Bcool — Bleak + Bfb — Bbo — Bgnpl (52) 

The efficiency of feedback /e is defined as the thermal 
energy gained (or lost) by the ISM (by the hot phase) at the 

end of the feedback process, divided by the energy injected 
by SNe. Cooling and leak-out are not directly associated 
with the action of SBs, while energy losses by blow-out and 
snowplow act in decreasing the thermal energy of the ISM 
by depleting the hot phase. We then define /e as: 

. -Bfb — -Bbo — -Bsnpl , , 

/e = (53) 

It is very important to notice that this quantity is not con- 
strained to be positive: in particular situations the depleting 
action of blow-outs and snowplows can overtake energy in- 
jection; in this case the net effect of SN explosions is a loss 
of thermal energy more than a gain. 

3.3 Metal flows 

For each generation of stars a fraction y of the restored mass 
is composed by new metals that are continually injected into 
the ISM. We call M-^ (where i=hot, cold, ★ or halo) the mass 
of metals in the various components, and Z\ = M-^ /M; their 
metallicities. In the instantaneous recycling approximation 
the system of equations for the metals is: 

' M^„t = -^hot(Meool + Msnpl + Meak + /boMnt) 
+ (1 - /bo)[-ZcoldA^evap + (-Zcold + 2/)Mrest] 
M,^ld = Zhot(Mcool + (17 /bo)M,„pl) -I- ZhaloMnf 

(Mevap + Msf) (54) 
Mf = Zcold(Msf -Mrest) 

A^halo = -.^haloMinf -I- Zhot[Mleak + /bo (Mint + M.npl)] 
+ /bo[-^cold Mevap -|- (.^cold + J/)Mrest] 

In this case mass is not conserved, the source term be- 
ing yAfrest. In these equations it is implicitly assumed that 
metals are efficiently mixed within each component. This 
assumption is reasonable for the hot and cold phases, but 
is clearly wrong for stars. In other words, Zi, is the average 
metallicity of stars but not that of the last generation, which 



contributes to enrichment. As a consequence, we use ^eoid 
for the metallicity of the newborn stars, while the actual 
distribution of metallicity of stars can easily be obtained by 
computing the evolution Zcoid-Aif* over time. Mixing of met- 
als within the halo is another delicate assumption; if gas is 
blown out in form of cold clouds then mixing may be ineffi- 
cient. A more refined modeling will be required in realistic 
cases. 

3.4 Parameters 

It is useful at this point to sum up the parameters introduced 
in the model. Some of them arc connected to the theory of 
stellar evolution or with the choice of the IMF; we do not 
regard them as free parameters of the model. In this pa- 
per we fix their values as follows: M*,Bn=120 Mq, /rest=0.2, 
tiife=2.7 X 10^ yr, y=0.04. 

The following quantities have been introduced in the 
various steps of the feedback model, and should be regarded 
as free parameters: ad, Pshape, Bsi, /evap, Tlivap, /cool. 
It is worth recalling that /evap and TLvap will be deter- 
mined in paper II, and that Tcvap plays no role if SBs arc 
pressure-confined. Other parameters are /bo, max, if blow-out 
is required to be more efficient than the simple geometri- 
cal model of Equation 31, and /pds, that regulates the in- 
jection of energy after PDS and is relevant only in case of 
PDS confinement. The parameters Tc and mi are presently 
kept constant (we have verified that the solutions do not 
change much for reasonable changes in these parameters). 
This parameter space is to be considered as minimal: all 
the computations presented above arc only useful ordor-of- 
magnitude estimates, so many quantities could in principle 
be fine-tuned (with the aid of new parameters) to reproduce, 
for instance, the results of detailed simulations. 

Finally, the total mass Mtot, the volume V (or equiv- 
alently the total density jOtot=A'/tot/V or surface density 
Stot=2.ffeffPtot), the infall time iinf, the vertical scale-height 
-ffeff, the velocity dispersion of clouds ctv and the geometry 
of leak-out d (fixed to 1 in the following) arc the parameters 
connected to the system in which feedback acts. These pa- 
rameters are obviously determined by the characteristics of 
the galaxy; in particular, the galaxy-halo system will never 
be a "closed box" as naively assumed here, so Mtot will not 
bo a constant; ptot will be fixed by the gravitational well of 
the dark-matter halo and the amount of angular momentum 
retained by the gas; Uv will be determined by gravitational 
infall, dissipation and re-acceleration by blasts; H^g will be 
determined by Cv and the surface density of the galaxy; <inf 
will be determined by the cooling and infall times of the 
dark-matter halo. To avoid the modeling of the galaxy at 
this stage, we consider for simplicity these parameters as in- 
dependent. The regions of these parameter space relevant to 
galaxies will be determined once their large-scale structure 
is fixed. 



4 RESULTS 

The system of Equations 43, 52 and 54 has been integrated 
with a standard Runge-Kutta algorithm with adaptive step- 
size (Press et al. 1992). The adaptive stepsize is computed 
considering only the equations for the mass and energy flows. 
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As mentioned above, the numerical integration is often del- 
icate, especially when the system switches from one regime 
to another. For this reason it is important to interpolate 
smoothly between different regimes. 

In the following wo fix the total mass of the system to 
10^^ M0. All the results can be simply rescaled to any total 
mass, as long as the mass allows for the presence of at least 
one collapsing cloud in the volume. The initial conditions 
are set by putting most gas into the halo, which is both 
physically reasonable and computationally convenient. The 
volume V is anyway filled with a small amount of mass in 
both cold and hot gas and a tiny amount of stars. This sets 
the system into a physically acceptable transient regime, al- 
lowing a smooth integration. We have verified that in general 
the precise choice of these initial conditions does not influ- 
ence the result as long as the system staxts in a way which 
is not pathological. Moreover, we set all primordial metallic- 
ities to 10~*. We specify the density of the system through 
the quantity ptot=Aftot/V^, i.e. the density that the system 
would have if all the mass were in the volume V and not in 
the halo. Clearly, the actual density of the ISM will always 
be smaller than ptot • For an easier comparison to astrophys- 
ical data, we show results in terms of the surface density 

Stot=2ptot^^eff- 



4.1 Feedback regimes for a reference choice of 
parameters 

We choose a reference set of parameters by fixing them to 
the typical (or tentative) values quoted above in the text: 

acl=2, ^61=1, /*=0.1, /cool=0.1, /Ushape=0.2, /evap=0.1, 

Tevap=10^ K. For this set of parameters we run the system 
of equations for a grid of values in the -ffeff-ptot (or equiv- 
alently -ffeff-Stot) plane, with a time scale of infall finf = 10* 
yr, which is suggested for the Milky Way by galaxy evo- 
lution models (see, e.g., Chiappini, Matteucci & Romano 
1997), and a velocity dispersion of clouds av — lO km s~^, 
typical of spiral discs. We stop the integration at three in- 
fall times, and check the actual regime of feedback. We have 
verified that the feedback regime may change with time, but 
in most CEises it does not change from 1 to Stini- 

Fig. 3a reports the regions in the _ffe£F-Stot plane in 
which various feedback regimes arc found. At low Etot and 
HeS values the SBs are able to blow out, and when they do 
they are always in the adiabatic stage. For increasing Hgs it 
is more and more difficult to blow out, so that SBs are kept 
confined by the external pressure in the adiabatic stage. The 
two regimes are roughly separated by the relation shown in 
the figure: 



Etot — 8 



1000 pc 



Mq pc" 



(55) 



At densities lower than this limit by roughly two orders 
of magnitude the system gets into a critical behaviour, where 
the hot phase is strongly depleted and the filling factor of the 
cold phase becomes high. This regime will be described in 
Section 4.3. At very high surface densities PDS confinement 
is met; this will be discussed later. 

If the vertical scale-length of a disc is set by dynamical 
equilibrium then Hgs = (7v/7rGEtot- Galaxy discs will then 
lie on the continuous lines shown in Fig. 3 for cTv=5, 10 and 



50 km s"^**. Discs with the canonical of a^—l km s~^ will 
be in the adiabatic blow-out regime, except possibly in the 
inner parts (where they blend with bulges) , while discs with 
(Tv>10 km s~^ will be in the adiabatic confinement regime. 
Bright spheroidal galaxies roughly follow a relation of the 
kind Re = 22(Mtot/10^^ M©)" '' kpc (proposed by Chiosi & 
Carraro 2002), that can be extrapolated to meet the locus 
of globular clusters. Identifying the effective radius Re with 
Hes, this relation is shown in the J?eff-Etot plane as a dashed 
lino. It is clear that feedback in a spheroid will typically be 
in the adiabatic confinement regime. 



4.2 Some examples 

We show here examples of the evolution of the system in 
various regimes. In particular we show the evolution, up to 
10 Gyr, of masses and metals of the four components, mass 
and energy flows, ISM and cloud properties. The evaporation 
and restoration rates, Mcvap and Mrcst , are not shown in the 
flgures as they are simply proportional to (and smaller than) 
the star- formation rate Mb!- 

Fig. 4 shows a Milky Way-like system in the adiabatic 
blow-out regime, with (ffeff, Ptot)=(100 pc, 0.1 Mq pc~^) 
or Etot=20 Mq pc"^ (it is denoted as MW in Fig. 3). Gas 
cools but does not transform promptly into stars; in fact, 
it is continually recycled into collapsing clouds. The final 
M^oid/Mt ratio is ~0.1. The hot phase regulates to a frac- 
tion Fh~2 X 10"''. The halo gas is continually recycled by in- 
fall, blow-out and leak-out. Regarding metals, the hot phase 
is promptly enriched to nearly solar values, followed by the 
halo gas, which receives the blown-out and leaked-out met- 
als; the final metallicity of the cold gas is solar, those of hot 
and halo components are 60 per cent higher. The star forma- 
tion rate Msf after a rise of ~1 Gyr decreases exponentially 
from ~20 to ~3 M© yr"'^ in ~9 Gyr; the average value of 
the star- formation rate is ~10 Mq yr~^. The infall rate Minf 
after ~2 Gyr regulates to slightly higher values, while the 
leak-out rate Mieak is very similar to the star-formation rate. 
Notably, leak-out dominates over blow-out, and the cooling 
term Mcooi is very low. The energy equation is characterized 
by the near equality of the two main flows, Efh and i^ieak; 
this reflects in the very stable value of Th~10'' K. The ISM 
is sclf-rcgulatcd and weakly varying over many infall times; 
its properties are nh~10~^ cm~^, nc~10 cm~^, P/fe~10^ K 
cm~^, and change by a factor ~5 from 1 to 10 Gyr. The 
porosity Qsb is low, indicating that active bubbles (i.e. ex- 
panding blasts) do not dominate the volume. However, if we 
assume that SBs are recognizable for a time tufc (see Section 
2.7), the porosity of "observed" SBs results as high as ~1 at 
a few Gyr and ~0.1 at the end of the integration; this implies 
an apparently bubble-dominated ISM. The flUing factor /c 
of the cold phase is ~0.1, while the fraction /coii of cold gas 
in collapsing clouds is slightly lower. The population of col- 
lapsing clouds is also rather stable, with masses in a range 
of roughly a factor of two around mcc~10^ Mq. The coagu- 
lation time tcoag is higher than the dynamical time tdynby a 



■* Here we have assumed that the vertieal seale-length of the 
molecular clouds is which amounts to halving H^g to take 

into account the easier blow-out of SBs that are off the mid-plane; 
see the discussion in Section 2.7 
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Figure 3. Feedback regimes in the i^cfl^^tot plane after three infall times. Filled circles denote adiabatic confinement, stars adiabatic 
blow-out, crosses PDS confinement. Critical cases are marked by circles. Thin lines separate regions with different regimes (see Equa- 
tion 55), thick continuous lines show the position of spiral discs with (starting from the lowest) (Tv=5, 10 or 50 km s~^, while the thick 
dashed line shows the average position of bulges. The labels MW and EL are relative to the examples shown below. Panel (a): reference 
choice of parameters. Panel (b): i?5i=0.3 and /cvap=0.05. 



factor of a few, and is ~5 x 10^ yr, not very different from 
the dynamical time of the Milky Way disc. Finally, both the 
number of clouds and the number of SNe per cloud are high 
enough to justify the assumptions of the model. 

Fig. 5 shows an example of adiabatic confinement (EL) 
that lies near the elliptical line (Fig. 3), with {HcS, ptot) = 
(3 kpc, O.3M0 pc"^) or Etot = 1800 M© pc~^. In this case, 
despite of the higher density, gas is consumed more slowly 
than the previous case, and the final fraction of AIcom/M^ 
is still ~0.5. The fraction of hot to cold gas is Fh~10~"^. 
The pattern of chemical enrichment is similar to the previ- 
ous cases, although metallicities are lower at the final time 
(due to the lower amount of gas consumed). Mass fiows peak 
to slightly lower values, and decrease less steeply at later 
times. Again, star formation regulates nearly to the infall 
rate after a few infall times, and leak-out is only slightly 
lower. Blow-out is obviously absent and cooling is again 
negligible. As before, feedback and leak-out energy fiows 
nearly compensate each other. The ISM is characterized by 
high pressure (P/fe~10* K cm~'^), density of both phases 
(nc~10^ cm~^, nh~10~'^ cm~'^), temperature of the hot 
phase (Th^l.S x 10^) and a correspondingly lower filling fac- 
tor of the cold phase (~2 x 10""^). Collapsing clouds are very 
small (mcc~3 x lO^M©), and this is due to the very high den- 
sity of the cold phase with the consequent low Jeans mass. 
Besides, the range of collapsing masses is tiny due to the low 
dynamical time and the consequent inefficient coagulation. 
This refiects into a low fraction of cold mass in collapsing 
clouds (/coii~10~^) and a low porosity of active SBs. 

This example is useful to understand the change in the 
behaviour of the system from the adiabatic blow-out to the 
adiabatic confinement regimes, but clearly a spheroid forms 



on shorter time-scales than 1 Gyr, and this leads (within 
the same physical time) to a consistently more rapid star 
formation, higher star- formation rates, higher enrichment, 
higher pressure and densities of ISM, slightly smaller collaps- 
ing clouds; however, the temperature of the hot phase and 
the filling factor and porosity of the cold phase are rather in- 
sensitive to the infall time. A similar trend is observed when 
the density is increased at fixed Hcs- 

The very low number of SNe per cloud in the adia- 
batic confinement case highlights a limit of applicability of 
the model in this case. However, it is easy to check that 
in the adiabatic confinement regime all terms in the sys- 
tem of Equations (43, 52 and 54) are independent of the 
actual size of SBs. The thermal energy of the first blast will 
be radiated away before the SNR manages to destroy the 
star-forming cloud, so a lower effective value of E51 will be 
reasonably used. The high value of the density of the cold 
phase highlights another problem. The reason why cold gas 
is not promptly consumed into stars is that it waits to be in- 
cluded in collapsing clouds. But for such high ric values the 
assumption that gravitational collapse is required to trigger 
the formation of H2 is probably wrong, and star formation 
is likely to be spread throughout the cold phase. This can be 
reproduced simply by forcing /coii to be unity; in this case 
the evolution of the system becomes trivial, the main mass 
fiows (A/sf, Minf, Mioak and Mevap) bccomc all proportional 
to each other and decay exponentially over one infall time. 

As shown in Appendix C (Fig. C2), to reach the PDS 
stage at Th=W^ K it is necessary to have rather high densi- 
ties Tih and relatively high mechanical luminosities Lag; the 
constraint tightens considerably at higher temperature. At 
high densities the Jeans mass is rather low, then Lag values 
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Figure 4. Evolution of the system for the MW adiabatic blow-out example with (-ffcffi Ptot)=(100 pc, 0.1 Mq pc~'^), Stot=2OM0 
pc~'^. The panels show the main properties of the systems, the quantities and their units are given in the labels. Time is linear, all the 
quantities given in ordinate are logarithmic. In the lower panels of mass and energy flows, the Af^f and Afi^ak curves and the Efb and 
^'loak curves are very similar and hardly distinguishable. 



are much smaller than unity, so that SBs are mostly kept 
confined in the adiabatic stage. Higher mechanical luminosi- 
ties could be achieved by increasing E^i , but this implies also 
a higher Th. As shown in Table CI, the ratio between tpds 
and tconi is proportional to L,^g ; as a consequence, 

the advantage in decreasing Lss is over-compensated by the 
increase in T^; as a consequence, PDS confinement is more 
easily achieved by lowering E51. We recall that for small 
clouds an effective lower E51 value is reasonably obtained 



because the thermal energy of the first blast is lost before 
the blast gets out of the collapsing cloud; moreover, for such 
dense clouds a lower value of /evap is likely (paper II). Fig. 3b 
shows the regimes in the i/eff-Stot plane for i?5i=0.3 and 
/cvap=0.05; while the limit between adiabatic blow-out and 
confinement (fig. 3a) is unchanged (and critical cases are 
found at densities higher by a factor 3), at densities roughly 
higher than 
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Figure 5. As in Fig. 4, for the EL adiabatic confinement example with {Hgff, ptot ) = (3 kpc, O.3M0 pc-3) or Stot=1800 M0 pc . 



Etot = 2500 (J/eff/1000 pc)-°-^ Mq pc"^ (56) 

PDS confinement is aciiieved. The reason wiiy it takes place 
preferentially at liigh H^g values is because leak-out is inef- 
ficient in depleting the hot phase and thus rih is higher. PDS 
confinement solutions are found also in Fig. 3a, though at 
very high surface densities. 

The evolution of the system depends sensitively on how 
the energy of SNe exploding after tun is given to the ISM. If 
mixing of hot phase and hot SB gas is slow, energy is pumped 
efficiently into the hot rarefied medium of the stalled bubble; 
this corresponds to /pds = l. Fig. 6 shows again the EL exam- 



ple of fig. 5 in this case. PDS confinement starts after a pe- 
riod of ~ 1.5 infall times of adiabatic confinement. Due to the 
sudden lower injection of energy into the ISM, Th decreases 
by nearly an order of magnitude at the start of PDS con- 
finement. This cooling has the effect of increasing tpds', the 
system then self-regulates to a configuration in which SBs 
stop just after PDS, so that the shell never acquires much 
mass. Pressure and densities are still high, but the filling fac- 
tor of the cold phase is as high as ~10~^. As a consequence 
collapsing clouds are bigger, and /coii is higher. At the on- 
set of PDS confinement the star-formation rate jumps to a 




Figure 6. As in Fig. 4, for the EL PDS confinement example (£51=0.3, /evap=0.05, /pds=l)- 



value of ~30 Mq yr~^, and then decreases slowly. Due to the 
lower Th, the fraction of hot gas is slightly high, fh~10~"^. 
In this regime cooling is more important than leak-out in 
terms of both mass and energy flows; we have verified that 
this is always the case for high Heg and Etot values. On the 
other hand, the snowplow flows are small (they are below 
the range of the energy flux panel) , indicating that the most 
relevant effect of PDS is on the structure of the ISM more 
than on the mass flows. Notably, the porosity of the SBs 
increases by more than an order of magnitude. At ~6 Gyr 



the ISM amounts roughly to 40 per cent of the total mass, 
and the solution switches back to adiabatic confinement. 

For smaller infall times or higher densities the PDS 
regime is triggered at an earlier time; again we have con- 
sistently higher star-formation rates, metallicities, pressure 
and densities of the ISM, while Th, /coii and Qsb are hardly 
affected. 

Fig. 7a shows the efficiencies of feedback /e (Equa- 
tion 53) in the examples discussed above. Unity values are 
obtained in the adiabatic confinement regime, while in the 
PDS confinement example the efficiency is just a few percent 
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Figure 7. (a) Efficiency of feedback /e as a function of time for the examples MW and EL with adiabatic or PDS (self-regulated or 
critical) confinement, (b) Star-formation time-scales (in units of the infall time) for the four cases of panel (a), (c) /e at three infall times 
(or at the onset of critical behaviour) for the cases shown in Fig. 3a (with the standard choice of parameters); triangles, crosses, stars, 
circles and asterisks are relative to _f/eft=10i 30, 100, 300 and 1000 pc; larger points denote the adiabatic blow-out regime, (d) Relation 
between star-formation time-scale and surface density; symbols as above, (e) Prediction of the Schmidt law versus the observational 
relation of Kennicut (1989); symbols as above. 



below unity (as apparent in panel (c) smaller values like 0.5 
are obtained at higher densities or lower infall times) . In the 
adiabatic blow-out regime most energy is lost to the halo 
and the efficiency is slightly lower than 0.1. Another way to 
quantify the efficiency of feedback is through its effect on 
star formation. We quantify it by the ratio between the star 
formation time t^f (Equation 38) and the infall time tinf. As 



shown in Fig. 7b, a higher efficiency corresponds to a longer 
star-formation time scale, but the correspondence is only 
qualitative. For instance, the EL PDS and adiabatic confine- 
ment cases have very similar feedback efficiencies but, due to 
the different ISM structure, star-formation time-scales that 
differ by more than a factor of two. 

It must be kept in mind that /e refers to the efficiency 
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with which the energy of SNe is given to the ISM. In the case 
of adiabatic blow-out, while ~5-10 of the energy is given to 
the ISM, ~5 per cent is lost in the destruction of the star- 
forming cloud (see paper II), and a comparable amount is 
likely lost in the acceleration of the bubble at blow-out. The 
remaining ~80 per cent of the budget will be available to 
heat up the halo gas. 

Fig. 7c and d show /e and tsf/imf at three infall times 
(3 Gyr, or at the final time in critical cases) for the grid of 
models shown in Fig. 3a. The efficiency of feedback /e jumps 
from a value 0.05-0.1 in the adiabatic blow-out CEises to 1 in 
the adiabatic confinement cases and then down to ~0.3 in 
the PDS confinement cases. The star-formation time-scale 
tsC is roughly fit by a relation: 



tgf — 25 £inf 



1 M0 pc" 



1 kpc 



(57) 



Critical and PDS confinement cases fall out of this relation. 

In spiral galaxies the star formation rate is well corre- 
lated with the amount of cold gas, following the so-called 
Schmidt (1959) law, quantified by Kcrmicut (1989) as Egf = 
(2.5 ± 0.7) X W-'^T.^-^^f^-^^ M© yr"^ kpc'^ (with Ecoid the 
surface density of cold gas in Mq pc~^). Fig. 7e shows the 
predictions of this relation for the same grid of models at 3 
Gyr, compared to the Kennicut relation. While the slope is 
accurately reproduced, the normalization depends on Hch, 
and is well reproduced for J7eff~50-100 pc. Bright spirals 
are known to have roughly constant surface densities and 
velocity dispersion of clouds, so the average value of Hcn is 
also constant and of order of its value in the solar neighbour- 
hood. So, the predictions of this model satisfy the Schmidt- 
Kennicut law also in its normalization. We have verified that 
this agreement holds in a very broad range of cases and at 
all times. This implies that this is a robust prediction of 
the model, but cannot be used to fine-tune the parameters. 
The Schmidt-Kennicut law is naturally obtained if star for- 
mation depends on the mass of the cold gas divided by its 
dynamical time. In our case (Equation 37) the relation is not 
built-in, due to the presence of the /coii fraction and to the 
fact that the dynamical time is computed on the actual and 
not average density of the cold phase. The Schmidt-like law 
then follows from the approximate constancy of /coii and fc- 

Finally, the Schmidt-Kennicut law is not followed in the 
external regions of spiral galaxies, where star formation is 
quenched. This is not predicted by the present model. How- 
ever, such star-formation edges are usually though to be an 
effect of differential rotation or, according to Schaye (2004), 
of photo-heating by the cosmological UV background. Both 
processes are not included here, so this disagreement is ex- 
pected and is not considered as a worry. 



4.3 Critical examples 

As long as the system of equations described in Section 3 
holds, the physical system is self-regulated; as we have seen, 
equilibrium solutions are found in which the ISM is rela- 
tively stable for many infall times. However, there are critical 
cases where the conditions for the existence of a two-phase 
medium are violated and the system of equations does not 
hold any more. 

At densities roughly two orders of magnitude lower than 



the limit shown by Equation 55 the system gets into a crit- 
ical regime, where the hot phase is strongly depleted and 
the filling factor of the cold phase fc is larger than that of 
the hot phase /h, violating the assumption of a cold phase 
fragmented into well-separated clouds. The cases where this 
happens are highlighted by a circle in Fig. 3. The reason for 
this behaviour is simple: at such low densities star formation 
is very weak, while for a thin structure leak-out is strong, so 
the hot phase cannot be sustained. This is not very informa- 
tive, as such low-density thin systems, if they exist, would 
be kept ionized by the cosmological UV background, so star 
formation would never start. 

A similar phenomenon happens in other cases for differ- 
ent reasonable choices of the parameters, with the difference 
that high values of fc are obtained not from the start but 
after some time. If we allow /bo, max to be as large as 0.8, 
the blow-out mass flux is much stronger and dominates over 
leak-out. For Etot~200-1000 M© pc~^, when Hes increases 
above ~1 kpc the system does not go into the adiabatic con- 
finement regime; SBs get bigger and blow-out gets stronger, 
severely depleting the hot phase. Pressure and densities are 
low, collapsing clouds are very massive (up to 10^ M©) and 
Qsb is high (even higher than unity). The hot phase is so 
strongly depleted that eventually, after roughly one infall 
time, /h</c- When the cold phase percolates the volume, the 
most likely outcome is collapse (the total mass of cold gas is 
surely larger than mj) and a sudden burst of star formation 
without any obvious external trigger. In other words, star 
formation would switch from a "candle" -like to a "bomb"- 
like solution. For a dynamical time tdyn~5 x lO'^ yr (in this 
case nc~l cm" "3) and for an efficiency /*~0.1, 10^° M© of 
cold gas would give rise to a starburst of tens of M© yr~^. 

Another example of critical behaviour is found when a 
high-density system goes into the PDS confinement regime 
(see the previous subsection). If the mixing of hot phase and 
hot internal gas is fast, then more thermal energy is lost to 
radiation and /pds=0. 737(1 - (ifin/tpds)'^'^) +0.263 (Equa- 
tion 51). In this case, when PDS confinement is triggered 
Th decreases dramatically, and the hot phase collapses in a 
very short time. Fig. 8 shows the evolution of the same ex- 
ample of Fig. 6 for this choice of /pds- The collapse of the 
hot phase takes place at the start of PDS, after 1.5 infall 
times (1.5 Gyr). The cold gas will be promptly consumed 
into stars; for a conservative value of /*=0.1 this will give 
rise to a brief star formation episode with Msf in excess of 
100 M© yr~^. The porosity of SBs takes values larger than 
unity at PDS, indicating the formation of a unique super-SB 
that will plausibly remove all the ISM not consumed by star 
formation. This super-wind will interact with the external 
halo gas, so that further infall will be halted for some time. 

The actuality of these critical solutions is uncertain, 

as they could be due to some of the simplifications intro- 
duced. In low-density cases, magnetic fields or turbulence 
could keep the cold medium fragmented even in the presence 
of low thermal pressure, or the hot phase could be replen- 
ished by mass fiows neglected here, while in high-density 
cases the existence of critical solutions depends sensitively 
on the way the energy of SNe exploding after the SB stalls 
is given to the ISM. In any case, the idea of a critical ISM 
deserves further investigations. 
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Figure 8. As in Fig. 6, but witli /pda=0. 737(1 - (tfln/tpds)"^'^) + 0.263. 



4.4 Probing the parameter space 

In the following we give a brief account of the effects of 
changing the various parameters within reasonable limits. 

The value of ctv enters formally in the coagulation time 
(Equation 11), so a hypothetical increase of ctv at fixed Hc« 
and Etot reflects in a decrease of the coagulation time, and 
an increase of the largest collapsing mass rriu (Equation 13) , 
with a corresponding increase of the collapsing mass rucc 
(Equation 14), the fraction of cold gas available for collapse 
/coll (Equation 15) and a resulting higher (and more quickly 
decreasing) star formation rate (Equation 37). However, the 



collapsing mass tticc is set mostly by the Jeans mass, so the 
results are rather insensitive to the precise value of CTv. Of 
course, in real systems a change in ctv would imply a change 
in //eft or Etot, with the known effects. 

Decreasing tinf results in a correspondingly stronger 
star-formation rate, and in a faster recycle of materials, 
while increasing it has an opposite effect. Again, the feed- 
back regimes do not change much; with fast infall adiabatic 
confinement is reached at slightly lower densities. 

The effect of decreasing a^i is that of moving mass to 
the high-mass end of the mass function of cloud, and thus 



Physical regimes for feedback 19 



to increase /coii and Msi- An increase of the star-formation 
rate is obtained also by increasing /* (Equation 37). How- 
ever, the two cases are rather different in terms of cold gas: 
with low /* and Oci values the cold gas is reprocessed by 
collapse many times, while with high values it is locked in 
the small clouds until it is processed by star formation. On 
the feedback regimes, lowering ccci or /* has the effect of 
moving the limit for adiabatic blow-out at lower densities, 
especially at low Hes values, and vice-versa for an increase 
of Oci or /* . 

An increase of /Jshape to unity lowers the Jeans mass 
(Equation 6), making blow-out more difficult; also, a lower 
mnnbor of SNc per cloud in the adiabatic confinement case 
is obtained. The opposite happens for a decrease of /ishape, 
which is justified if collapsing clouds are supported by ki- 
netic or magnetic pressure; for //shape=0.05 adiabatic blow- 
out is easily reached at densities higher by and order of 
magnitude than the limit given in Equation 55. 

As mentioned above, a low value of E51 leads to lower 
Th, lower pressure and a corresponding increase of the Jeans 
mass; moreover, PDS confinement is found at high H^g and 
Stot values. Increasing E51 leads to an increase of Th, that 
for i?5i=10 can reach extremely high and unlikely values. 
The limit between the adiabatic blow-out and confinement 
regimes does not depend much on £51. 

As cooling is a relatively modest mass fiux, the results 
do not depend sensitively on the parameter /cooi, with the 
exception that PDS confinement is reached more easily for 
E5i=l whenever /cooi<Cl. 

The parameters /evap and Tevap mostly affect the adiar 
batic blow-out regime; in particular Tevap influences strongly 
the resulting Th. A low value of /evap increases the number 
of critical blow-out cases, until Tevap is increased to compen- 
sate for the lack of evaporated mass. If both parameters are 
increased, the limit for adiabatic blow-out lowers. Finally, 
with a high /evap value PDS confinement is reached at high 
densities even for £51 = 1. 

4.5 On the vertical scale-height of the hot phcise 

The assumption of one single vertical scale-height for both 
cold and hot phases is clearly over-simplistic; it can be re- 
laxed by assuming two different scale-heights, but in absence 
of further constraints this new degree of freedom would not 
contribute significantly to the understanding of the problem. 
Anyway, the assumption is sensible as the hot gas is contin- 
ually replenished within i?eff- The hot gas that leaks out 
into the halo is likely to settle in a low-density layer that 
surrounds the galaxy. Such layer is observed in the Milky 
Way (see, e.g., Jenkins 2002) as well as in nearby galaxies 
(see, e.g., Ferguson ct al. 1995; Fraternali et al. 2002). The 
presence of a sufficiently steep decreasing density gradient at 
Hes is enough to guarantee that this layer does not hamper 
the blow-out of SBs. Indeed, the interaction of blown-out gas 
and such a layer could be at the origin of the observed cor- 
relation of X-ray and Ha fluxes in nearby starburst galaxies 
(Strickland et al. 2002). 

We can estimate the thickness of this layer, that we call 
Hh, as follows. We assume that the expansion of the leaked- 
out gas stops after one cooling time. For a gas in adiabatic 
expansion we have n'-^ = nh-ffetr/^^h, Th = Th(i?eff /^fh)^~^ 
and Cg h = Ca,h{Heff/Hh)^''~^^^^, where the prime indicates 



the quantities relative to the gas layer and 7 is the adiabatic 
index, assumed to be 5/3. It is easy to see that, for the 
cooling function given in Equation 23 the cooling time is 
constant with Hh- The equilibrium value of Hh will satisfy 
the condition fcooi = Hh/c[. To compute it, we take into 
account that there is continuous injection of hot gas with a 
roughly constant rate within one cooling time. We obtain: 

//h = ^i/eff(^2^)'^'. (58) 

For the adiabatic blow-out case of Fig. 4 we obtain Hh~2 
kpc. This will be an overestimate, as gravity is likely to be 
important at such distances from the centre of the galaxy. 

Alternatively the expansion of the hot phase could 
progress in such a way not to create a steep density gradient. 
In this case the hot phase would be contained in a layer of 
thickness H^ and nearly constant density. Exploiting Equa- 
tion 58, it is possible to include the dynamical evolution of 
such a layer in the system of equations; in this case leak-out 
would not be considered as the hot phase formally never gets 
back to the halo. As a result, the hot gas produced by SNe 
pushes the layer to high enough Hh values that adiabatic 
confinement is always achieved. The predictions of this ver- 
sion of the model are found in striking disagreement with 
observations in the Milky Way-like case (see next section), 
so the possibility of a hot gas layer with a roughly constant 
density will not be further considered. 



5 DISCUSSION 
5.1 The Milky Way 

To highlight the predictive power of the model, we compare 
the results of the Milky Way-like adiabatic blow-out case 
shown in Fig. 4 with available observational evidence. It is 
useful to stress that no accurate modeling or fine-tuning 
of parameters is attempted at this stage, so an order-of- 
magnitude agreement is considered a success. The predic- 
tions of rih (~10-^ cm-3), ric (~10 cm'^) and Th (~ 3 x 10'^ 
K) are in broad agreement with the ISM of the Milky Way. 
Thermal pressure (-10^ K cm •') is in line with obser- 
vations, but is an order of magnitude lower than the ob- 
served total pressure, which is dominated by turbulent and 
magnetic contributions. The mass ratio of cold gas to stars 
(~0.1) is correctly reproduced after 10 Gyr. 

The star formation rate Msi is slowly decreasing in time, 
with a ratio of average to final rates of ~3. This roughly 
consistent with the results of the chemical evolution model 
of Chiappini et al. (1997). 

Both the average value and the range of the masses of 
collapsing clouds are smaller than those observed for molec- 
ular clouds, that can be as large as md > 10^ M©. Large 
collapsing clouds are easily obtained by using a very low 
value for /ishapo, on the ground that kinetic support deter- 
mines the Jeans mass mj (Section 2.3). The small range of 
collapsing cloud masses reflects in the low values of /coii , the 
fraction of cold mass in collapsing clouds; it is predicted to 
be ~5 per cent, at variance to the observed ~50 per cent. To 
obtain higher /coii values it is useful to decrease ad to the ob- 
served value of 1.6, decreasing also /« to avoid excessive star 
formation (this is also consistent with observations). How- 
ever, good /coll values are obtained only by allowing clouds 
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to coagulate for at least 10 dynamical times. This unreal- 
istic value is not worrisome if we consider the uncertainty 
connected to the coagulation picture; while the coagulation 
time, icoag~5 x lO'^ yr, is coincident ally similar to the time 
interval between the sweeping of two spiral arms, the for- 
mation of molecular clouds in the converging flows of spiral 
arms (Ballesteros-Paredes, Vazquez-Semadeni & Scalo 1999) 
could easily be more efficient than random aggregations of 
clouds, and this could be mimicked by allowing coagulation 
to work for many dynamical times. 

The thickness of the layer of leaked-out hot gas is pre- 
dicted to be ~2 kpc, in rough agreement with the value of ~3 
kpc estimated for the Milky way by Savage ct al. (2000; sec 
also Jenkins 2002). However, this estimate is based on FUSE 
detection of O VI absorption lines of OB stars; this method is 
sensitive to temperatures in a narrow range around ~3 x 10^ 
K. The gas leaking out at Th=10® K has adiabatically cooled 
to ~1.5 X 10^ K at 2 kpc, while the temperature of ~3 x 10^ 
K is reached at ~0.6 kpc, significantly less than observed. 
However, this prediction depends sensitively on the param- 
eter Tevap, that influences the temperature of the hot phase. 
If this parameter is increased by a factor of 3 (a reasonable 
choice according to paper II), the resulting layer is predicted 
to be 10 kpc thick, reaching a temperature of ~3 x 10^ K at 
~3 kpc as observed. 

With reasonable choices of the parameters, and allow- 
ing coagulation to work for 10 dynamical times, it is possible 
to reproduce all these properties of the Milky Way. The rea- 
son why we do not stress this result is because we consider 
the present model too simple to draw significant conclusions 
from it. By interfacing this model with an algorithm for disc 
formation in a cosmological dark-matter halo and including 
the effect of differential rotation and spiral arms it will be 
possible to produce accurate predictions for the Milky Way, 
including galactic fountains, high-velocity clouds, chemical 
enrichment of the various components, chemical gradients 
along the disc and so on. By reproducing the observed Milky 
Way it will be possible to constrain most model parameters 
by modeling just one galaxy. 

5.2 Critical solutions and the triggering of 
galactic winds 

Although a proper modeling of galactic winds requires spec- 
ifying the properties of the dark-matter halo hosting the 
galaxy, it is interesting to analyse the cases in which feed- 
back could lead to the removal of a significant quantity of 
ISM from a galaxy. As removal of gas from a halo with low 
circular velocity can be achieved even with a single SB (see, 
e.g., Fcrrara & Tolstoy 2000), we will concentrate on bright 
galaxies. Blow-out leads to the expulsion of matter with a 
velocity that, in the example of Fig. 4, ranges from ~250 
km s~^ at 1 Gyr to more than 500 km s~^ at later times, 
so if these clouds are not slowed down significantly by the 
halo gas (e.g. by the layer of leaked-out hot gas) they may es- 
cape even from relatively high-mass halos. Anyway, blow-out 
flows are never very strong, so blow-out is unlikely to lead to 
massive removal of ISM from a galaxy. Besides, leaked-out 
gas cools below lO' K at ~10 kpc, so it will bo emitted as a 
tenuous wind from the low-mass halos, but will be retained 
by the halo of a bright galaxy. Leaked-out gas is much hotter 
in the adiabatic confinement regime, so this gas will be able 



to escape from moderate-sized halos, but will be retained for 
instance in big elliptical galaxies. In conclusion, as long as 
blasts propagate into the hot phase and SBs do not perco- 
late the volume, the removal of mass is inefficient in bright 

galaxies. 

This conclusion changes in the critical cases, where a 
significant amount of gas accumulated for some time is con- 
sumed in a few dynamical times, and when the porosity of 
SBs gets unity value. 

Critical solutions are found at least in three cases: (i) 
thin, very low-Stot systems in the adiabatic blow-out regime, 
(ii) thick, moderate- Etot systems in the adiabatic blow-out 
regime (in case of very efficient blow-out) and (iii) thick, 
high-Stot systems in the PDS confinement regime (with low 
E51 and low /pds). While case (i) has no astrophysical rel- 
evance (such systems would anyway be kept ionized by the 
cosmological UV background), case (ii) may correspond to 
some gas-rich dwarf galaxies and case (iii) to high-redshift 
spheroids. 

When these systems become critical, all the cold phase 
collapses and gives rise to diffuse star formation. For a con- 
servative /♦ value of 0.1, we estimated star-formation rates 
of tens of Mq yr""*^ for case (ii) and in excess of 100 M0 
yr~^ for case (iii). In such big bursts /* could well take unit 
values, thus boosting star formation rates even higher; on 
the other hand, if the transition from the "candle"- to the 
"bomb"-like regime is not as quick as assumed, star forma- 
tion rates will be lower. Analogously to what happens in 
star-forming clouds, the exploding SNe will propagate into 
the diffuse cold phase, going soon in the PDS stage and 
then percolating the volume. This will give raise to a unique 
SB with a very high mechanical luminosity, able to sweep 
the whole galaxy. This snowplow will eventually blow out 
of the galaxy and then fragment because of Raileigh- Taylor 
instabilities. If the momentum of the gas in the fragmented 
snowplow at this point is sufficient, it will be thrown out of 
the galaxy. 

Percolation of SBs gives a similar effect if it takes place 
in the PDS confinement (or blow-out) regime; though it has 
been assumed for simplicity that clouds pierce the snowplow, 
this is likely true only for the largest and densest clouds, 
that would likely be star-forming in this case. The effect of 
a percolation of collapsed shells would be to create a super- 
SB that sweeps the ISM, pushing part of the gas out of the 
galaxy in form of cold clouds while the rest is compressed 
toward the centre of the galaxy. This is found, for instance, 
in the simulations of primordial galaxies by Mori, Ferrara & 
Madau (2002); as the physics is the same, their conclusion 
can be extended to larger, lower-density galaxies, as long 
as percolation of SBs in the PDS confinement regime is ob- 
tained. Obviously, the gas concoutratcd at the centre would 
give rise to a secondary burst of star formation that would 
pump further energy into the super-SB. 

Percolation of SBs in the adiabatic stage is likely to 
have a smaller effect, as the blast would continue to prop- 
agate into the hot phase. The cold phase would be affected 
by the relatively inefficient processes of thermo-evaporation 
and cloud dragging (preferentially in the radial direction). 
These same processes arc in place also in presence of a hot 
phase that continually leaks out of the galaxy. In the adi- 
abatic confinement regime the cold phase is so dense and 
with such a low filling factor that these effect are likely to 
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be inefficient, while in the case of adiabatic blow-out from 
a disc the dimensionality of the system would presumably 
lead to a funneling of the energy in the vertical direction, 
thus making a massive removal of gas unlikely. 

As a matter of fact, the condition Qsb>l is never met 
in the adiabatic confinement regime, while it is achieved in 
the critical cases discussed above (Fig. 8). 

There are other cases in which the system may become 
critical. In presence of very high pressure the density of the 
cold phase can be so high (say >10'^ cm^^) that the forma- 
tion of H2 is triggered even in absence of collapse. If this 
limit is reached when much gas is accumulated, this may 
give raise to a sudden burst of star formation. 

Finally, a critical behaviour of the system can be trig- 
gered from outside. For instance, a strong tidal perturbation 
(or a merger) would act in two important ways, by lower- 
ing the Jeans mass of the clouds (because of the external 
pressure) and by thickening the structure, thus allowing a 
disc-like system to switch from adiabatic blow-out to adi- 
abatic confinement; this would lower the Jeans mass even 
more and decrease the dynamical time. The effect would be 
a rapid consumption of the accumulated cold phase and a 
likely percolation of SBs. 

5.3 Simplified models 

The present model can be generalized to reproduce the com- 
ponents of real galaxies, like disc, bulge and halo, and then 
interfaced with a galaxy formation code that includes the 
meiss assembly of dark-matter halos, cooling inside those ha- 
los, disc formation, galaxy mergers, interaction with galaxy 
clusters etc.. However, it is much more convenient to de- 
vice a set of approximate analytic solutions to this feedback 
model. These solutions can also be adapted to model the 
"sub-grid" physics of feedback in N-body simulations. 

The solutions in the ffcff-Stot plane can be divided into 
four main regions where different regimes are met (adiabatic 
blow-out, adiabatic confinement, PDS confinement, critical 
blow-out cases). These regions are separated by limiting re- 
lations of the kind Etot = Stot,o(^^off /I kpc)^""™, where the 
exponent is usually in the range 0.5 1. At the lowest 
densities (a factor 10^ lower than Equation 55 for the ref- 
erence choice of parameters) systems are critical, but they 
will most likely be kept ionized by the cosmological UV back- 
ground, so they will simply not evolve. At densities below 
Equation 55 (again for the reference choice of parameters) 
the system is in the adiabatic blow-out regime. Above that 
limit it gets into the adiabatic confinement regime. PDS 
confinement is reached at densities higher than Equation 56 
(valid for £51=0.8 and /cvap=0.05). 

For each regime a simplified solution can be obtained 
by noticing the following facts that are found to hold in 
most cases: 1) Th is nearly constant, and equal to a value 
that mostly depends on the regime and on E51; 2) Fh, /coii 
and /c are nearly constant to values that mostly depend on 
the regime; 3) cooling is negligible in all cases but those in 
PDS confinement, where leak-out is negligible; 4) leak-out 
dominates over blow-out in non-critical cases if /bo is not 
largo. With these assumption it is relatively easy to solve 
the system of Equations 43 for the mass flows, while typical 
values of Th, Fh, /coii and fc the different feedback regimes 
have been given in Section 4. However, a proper presenta- 



tion of these simplified solutions requires some discussion 
that is out of place in this paper, so they will be presented 
elsewhere. 

5.4 Limitations and further work 

The merit of such relatively simple modeling is to high- 
light the possible physical regimes one should expect once a 
more complete calculation is performed. However, there are 
a number of limitations that have to be carefully taken into 
account to assess the validity of the results presented hero. 

(i) The Sedov solution for the SBs in the adiabatic stage 
is only a rough approximation of reality. There is a long list 
of effects, mentioned above and in part described by Ostriker 
and McKee (1988), that influence the dynamics of SBs. How- 
ever, as long as the SBs expand in the relatively smoother 
hot component, it is likely that the Sedov solution gives the 
roughly correct evolution and functional dependences for the 
SBs. 

(ii) As already mentioned, thermal conduction at the in- 
terface of cold and hot phases can make part of the cold gas 
evaporate. It has been verified that the impact of thermo- 
evaporation is small in the mass flows even if it is not 
quenched by magnetic fields, as the thermo-evaporated mass 
is usually smaller than the evaporated mass of the collapsing 
cloud whenever /cvap is not much smaller than one. 

(iii) Type la SNe have not been considered. However, 
their introduction is straightforward in this model; they will 
interact directly with the ISM through a set of uncorrelated 
SNRs. As shown by Recchi et al. (2002), Type la SNe may 
be very important because they explode after a burst of star 
formation, and can contribute to maintain the hot phase 
when most cold gas is consumed. 

(iv) Cosmic rays are known to be in rough equipartition 
with turbulence and magnetic fields. They are accelerated 
by the shocks generated by SNRs and SBs, directly or in- 
directly through turbulence (see, e.g., Longair, 1981). The 
role of cosmic rays, which are confined within the galaxy by 
magnetic fields, is that they distribute their energy to all 
the ISM, and not only to the densest collapsing clouds. So, 
they could give an important contribution to the mass and 
energy flows. 

(v) There are other channels of mass and energy ex- 
change between components that we are not considering 
here. One is the decay of turbulence driven by the kinetic en- 
ergy of SNe, the other is the presence of a significant amount 
of mass in a warm phase, that can receive part of the en- 
ergy of the blast and radiate it. Also, UV light coming from 
massive stars or from an external UV background could be 
responsible (together with thermo-evaporation) for continu- 
ous evaporation of the cold phase. Although some analytical 
estimates of these effects are possible, accurate numerical 
simulations will be necessary to assess the importance of 
these processes. 

(vi) In realistic situations the ISM is subject to many 
influences, like spiral arms, differential rotation, tidal distur- 
bances, mergers, ram pressure from hot halo gas (like in el- 
lipticals or clusters) etc.. All these processes can be modeled 
once the global structure of the galaxy and its environment 
are specified. For instance, the passage of a spiral arm can 
be modeled by a periodic decrease of the Jeans mass due 
to the external pressure term. As in the adiabatic blow-out 
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case the coagulation time is tcoag~5 x 10*^ yr, of the same 
order of the frequency of spiral arms, the clouds have just 
time to coagulate to increase their mass by a factor of a few 
before the spiral arm sweeps again. So, a moderate decrease 
of the Jeans mass would suffice in guaranteeing that star 
formation takes place mainly in the spiral arms, even in ab- 
sence of a more explicit connection, like that proposed by 
Ballesteros-Paxedes et al. (1999). 

(vii) The model presented here is assumed to be valid 
in the regime where many generations of collapsing clouds 
self-regulate in forming a galaxy. For dwarf galaxies some 
changes in the model are necessary: firstly it is important 
to check that at least one collapsing cloud is present; sec- 
ondly, the first episode of star formation could itself cause a 
complete blow-away of the ISM (see, e.g., Fcrrara & Tolstoy 
2000) , so that the system may never get into a self-regulated 
regime. Moreover, for dark matter halos with small circular 
velocities the leaked-out or blown-out gas will most likely be 
lost to the inter-galactic medium. 

(viii) This model gives satisfactory predictions for the 
state of the ISM in a Milky Way-like situation. The high- 
density, bursting cases are subject to much weaker observa- 
tional constraints, not only for the paucity of very nearby 
starbursts but also for the presence of dust that hampers 
observations in the optical, UV and soft-X bands. Besides, 
the extrapolation of the assumptions that are successful at 
low densities is not straightforward. For instance, large col- 
lapsing clouds could be generated in cooling flows or, as 
mentioned above, star formation could be triggered even in 
clouds smaller than the Jeans mass when thermal pressure 
maJtes them denser than a threshold density at which H2 
starts to form. Careful comparison with available observa- 
tions is needed to constrain the parameters of feedback in 
the starburst cases. 

(ix) This model can be improved to give more accu- 
rate predictions on obscrvablcs related to the ISM. This 
would require proper (numerical) modeling of the intermedi- 
ate warm phase(s) and of the ionization equilibrium between 
the phases. 

5.5 Comparison with previous works 

The model of the ISM of McKee & Ostriker (1977) has been 
a reference model for years, although the picture based on 
compressible supersonic MHD turbulence is now emerging. 
The model presented here has many points in common with 
McKee & Ostriker (1977), but presents many improvements: 
(i) we address the dynamics of the ISM, including star for- 
mation and feedback; (ii) wo assume no equilibrium, but 
investigate on the conditions that lead to self-regulated or 
critical ISM; (iii) we take into account the correlated na- 
ture of Type II SNe. Besides, we do not consider the warm 
phase and its ionization equilibrium with the hot phase. In 
the present model we do not require unit porosity of SBs 
to justify its presence of a hot phase (which sometimes can- 
not even be maintained). As a matter of fact, unit porosity 
of SBs is hardly reached in non-critical cases; however, un- 
correlated adiabatic SNRs that stop at Rho' would have a 
porosity of order one in the MW example, but not in the EL 
one in the adiabatic confinement case. So, on the light of the 
present results, Q = 1 is at best an unnecessary assumption. 
Some recent works on galaxy formation by Silk (1997; 



2000), Efstathiou (2000), Shu, Mo & Mao (2003) or Springer 
& Hernquist (2003) present models of feedback and star for- 
mation based at least in part on the McKee & Ostriker 
(1977) model. In these cases Type II SNe are assumed to 
be uncorrclated and the ISM is assumed to be self-regulated 
to a unit value of the porosity of SNRs. For instance, Ef- 
stathiou (2000) fixes the star formation by assuming equi- 
librium between the kinetic energy acquired by cold clouds 
at shocks and that lost by coagulation, while Silk (1997, 

2000) connects star formation to the dynamics of the disc 
by requiring the Toomre Q-parameter to be unity and postu- 
lating the identity of the time scale for star formation with 
the viscous time scale: this is an important ingredient for 
obtaining exponential discs, but the nature of this identity 
remains unexplained. The model presented here does not 
assume any equilibrium, and does not use any ingredient of 
disc dynamics, thus being applicable in virtually all situa- 
tions. While it is clear that disc dynamics will influence the 
evolution of a spiral galaxy, our results suggest that most 
properties of galaxy formation can be understood simply as 
a chain of local processes. 

An alternative to the present modeling is to consider 
the ISM as turbulent. As shown by Avila-Reese & Vazquez- 
Semadeni (2001) the ISM can be considered as a globally 
turbulent medium, with turbulence forced in specific places 
(the star-forming regions) and propagating throughout the 
volume. The "diffusion" velocity of turbulence is connected 
to the time scale of decay of turbulence. Both this group, 
that uses a 2D code, and Mac Low et al. (1998), who use a 
3D code, find that turbulence decays as t~°' with a ~ 0.8. 
It is easy to show that (Avila-Reese & Vazquez-Semadeni 

2001) , when turbulence is forced in some specific sites, the 
nils velocity of turbulence scales with distance from the forc- 
ing region as Urma oc while the decay distance of 
turbulence grows with time as I oc t^'"^^. It is remarkable 
that for a = 0.8 the two exponents are exactly equal to our 
relations Rsh oc t"''' and Vsh oc R^" '^. While a direct phys- 
ical interpretation of this fact may be misleading without 
further investigation, it is clear that the propagation of en- 
ergy through the ISM by isolated spherical blasts is not in 
clear contradiction with the results of the turbulent model. 
This confirms the validity of a simple treatment as a first 
approximation. 



6 CONCLUSIONS 

We have presented a model for feedback in galaxy formation, 
based on a two-phase ISM, that does not restrict to self- 
regulated, equilibrium solutions and neglects (for simplicity) 
the global structure of the galaxy, apart from its density, ver- 
tical scale-height and velocity dispersion of clouds. From the 
dynamics of the SBs that arise from the collapsing "molecu- 
lar" clouds, we have identified four possible regimes of feed- 
back, depending on whether SBs blow out of the "disc" or 
remain pressure-confined, and whether they have time to en- 
ter the PDS stage. For a reference set of parameter values we 
have studied the dynamics of the system in the vertical scale- 
height surface density plane, identifying the regions of the 
plane corresponding to different regimes. Both blow-out and 
confinement mostly take place in the adiabatic regime. In a 
Milky Way-like adiabatic blow-out case, the main charac- 
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teristics of the ISM of the Galaxy arc broadly reproduced. 
In the adiabatic confinement regime the ISM is predicted to 
have higher pressure, temperature of hot phase and densi- 
ties of both phcises, and smaller collapsing clouds; in some 
cases the density of the cold phase could be high enough 
to trigger diffuse star formation. PDS confinement is found 
for high-density, thick structures in significant regions of the 
parameter space. In this case feedback is less effective, the 
hot phase cooler and star formation quicker. 

In many cases the system becomes critical, in the sense 
that the hot phase is severely depleted and the cold phase 
percolates the whole volume. This happens for very low- 
density thin systems (that would however be kept ionized 
by the cosmological UV background), in some regions of 
the parameter space also for low-density thick systems in 
adiabatic blow-out (that may correspond to some gas-rich 
dwarf galaxies) and for high-density thick systems in PDS 
confinement (that may correspond to high-rcdshift galaxies). 
The most likely result of this critical behaviour is the sudden 
consumption by star formation of the cold gas accumulated 
by the galaxy; the dynamics switches from a "candle"- to a 
"bomb"-like solution. 

The porosity of SBs is usually found to be much lower 
than unity. However, in some cases unit porosity is found 
while SBs are in the PDS stage. This corresponds to the 
formation of a super-SB that sweeps the whole galaxy, re- 
moving most ISM from it. These events, together with the 
critical solutions, are likely connected to the triggering of 
galactic winds. 

With respect to previous models of feedback, the main 
parameters that are typically present, as the efficiency of 
feedback, the Schmidt law with its normalization, or the 
rate of blow-out and leak-out of gas from a star-forming 
galaxy, are predictions of the present model. The parameter 
space is connected to the properties of the ISM, and can 
thus be constrained by observations of the Milky Way and 
nearby galaxies; most parameters can be fixed in principle by 
reproducing only the Milky Way. Moreover, the mass flows 
used in this model can be fine tuned by comparing with 
future detailed simulations of the ISM in a forming galaxy 
that include all the main physical processes though to be at 
work. 

This model does not restrict to self-regulated ISM, and 
presents a rich variety of solutions with a relatively limited 
set of parameters. Although the turbulent nature of the ISM 
is not explicitly taken into account, the model is thought to 
give a good approximation to the solution of the feedback 
problem. The feedback regimes found here can be used, to- 
gether with the refinements of the model that will be given 
in upcoming papers, to construct a realistic grid of feed- 
back solutions to be used in galaxy formation codes, either 
semi-analytic or numeric. 
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APPENDIX A: LIST OF FREQUENTLY USED 
SYMBOLS 



Oci Slope of mass function of clouds 

Oci Radius of cloud 

Cs,h Sound speed of the hot phase 

-E51 Energy released by a SN 

e'^^ Thermal energy of a SB 

^(kin) Kinetic energy of a SB 

Es-n Rate of energy release from SNe 

-Bcooi Rate of energy loss by cooling 

i^snpi Rate of energy loss by snowplows 

Eho Rate of energy loss by blow-out 

Sioak Rate of energy loss by leak-out 

Eih Rate of energy gain from SBs 

Eyioi Net energy flux of the hot phase 

/e Efficiency of feedback 

/h,/c Filling factors of the two phases 

/coll Fraction of cold gas in collapsing clouds 

fi, Efficiency of star formation 

/evap Evaporated fraction of collapsed cloud 

/bo Fraction of swept gas blown out by a SB 

/bo.max Largest value of /bo 

/cool Fraction of cooled gas in a cooling flow 

/rest Fraction of restored mass 

/pds Release of energy after PDS confinement 

Fh Fraction of hot gas 

Hed Vertical scale-height of the system 

Hes,h Dynamical vertical scale-height of hot g£is 

Hh Height of the layer of hot leaked-out gas 

L38 Mechanical luminosity of a SB 

mci Mass of clouds 

mi Lower cutoff mass of clouds 

ruu Upper cutoff mass of clouds 

mj Jeans mass of clouds 

mcc Mass of the collapsing cloud 

Msw Mass swept by a SB 

Mint Swept mass that is still hot inside a SB 

Mt,Bn Mass of formed stars per SN 

Mtot Total mass of the system 

Mi Mass of the i component* 

Mass of metals in the i component* 

Mi Net mass flux of the i component* 

M^ Net metal flux of the i component* 
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Mint 


Infall rate 


Mcool 


Cooling rate 


Mleak 


Leak-out rate 


Msf 


Star formation rate 


Mrest 


Restoration rate 




Evaporation rate 


Mint 


Swooping rate minus snowplow rate 


Msnpl 


Snowplow rate 




Molecular weights of the two phases 




Shape parameter for collapsing clouds 


TtVi .?7.r 


Density of the two phases 


N^i 


Mass function of clouds 




Total number of collapsing clouds 




Number of SNe in a collapsing cloud 


P 


Pressure of the ISM 


Qsh 


Porosity of SBs 




Radius of SB 


Rpds 


Radius of shell collapse for a SB 




Confinement radius of a SB 


Rho' 


Final radius of a blown-out unconfined SB 


Rfin 


Final radius of a SB 


Rsn 


Rate of SN explosions in a collapsing cloud 


Ptot 


Total density of the system 




Average densities of the two phases 


Cv 


Velocity dispersion of clouds 


S^ot 


Total surface density of the system 


tcoag 


Coagulation time 


dyn 


Dynamical time of clouds 


tcool 


Cooling time of hot gas 




Sound crossing-time of a SB 


tpdT 


Time of shell collapse for a SB 




Confinement time for a SB 


thfe 


Lifetime of an 8 M© star 


iho 


Time of first blow-out of a SB 




Final time of a blown-out unconfined SB 


ifin 


Final time of a SB 


t'lnf 


TTifall l^iTTip 'ipalp 


ileak 


Leak-out time scale 




Star-formation time scale 




Temperature of the two phases 


^evap 


Temperature of evaporated mass 


Vsb 


Velocity of SB 


2/ 


Yield from massive stars 


Zi 


Mctallicity of the i component* 


Cm 


Metallicity in solar units 


Notes: 





i = hot, cold, *, halo 



APPENDIX B: TIME-SCALES FOR THE 
COAGULATION OF COLD CLOUDS 

We obtain here the time-scales for coagulation, given in Sec- 
tion 2.4, from the Smoluchowski equation of kinetic aggre- 
gations. This demonstration follows Cavaliere et al. (1992). 
The Smoluchowski evolution equation for the mass function 
of clouds n{m; t) is: 



dn 
dt 



dm' K{m' , m — m')n{m'\ t)n{m — m; t) 



n{m;t) / dm,'K{m,m')n{m';t) . 



(Bl) 



The kernel for aggregations of clouds 1 and 2, K(rni,Tn2) 
is given by Equation 9, where the cross-section for inter- 
actions is given by Equation 10. We define a typical mass 
for the mass function (which is then identified with the 
upper cutoff ruu) and scale all masses with rre* through the 
unidimensional variable i' — m/m^. We separate the kernel 
for aggregation of clouds into geometric and resonant terms 
(the latter term being considered only in this appendix) , and 
write the two terms as follows: 



I 1/3n2\ 
geom "f-* W^l ^" ^2 ) )^ 



(B2) 



1/3 



1/3n 



We call A the exponent of m* in the equations. The two T 

functions are respectively J-^som = PcTt 

and .7^-63 = Pc27rG'(47rpc/3)~"'''^((vap)~^)v. Let's assume 

that the mass function is expressible as: 



n(m; t)dm = m^^{v)di> . 



(B3) 



This is valid if the slope of the mass function is fixed and 
if mi <C m-a. Inserting this anstatz into the Smoluchowski 
equation, and considering only the time-dependent terms, 
we easily obtain the equation: 



(B4) 



This is valid for the two coagulation modes. This equation 
admits the solution: 



m*(t) = m*o(l - (A - l)jrm^o"'(i " *o))'^''"^^ ■ 



(B5) 



We can then define a coagulation time as tcoag = 1 l^fn^o ^ ■ 
For the two coagulation modes we obtain: 



m*(t) = m*o ( 1 + 



t-to 

Sicoae 



m*(t) = m*o ( 1 



t-to 

3icoae: 



(GEOM.) 



(RES.) . 



(B6) 



It is easy to verify that the coagulation time is given by 
Equation 11 in the case of geometrical interactions, while in 

the resonant case: 



tc 



/47r\i/3 J___ 



2/3 Pc 1/3 



(B7) 



Finally, it is worth noticing that the slope of the mass 
function is assumed not to change during the evolution of the 
system, while coagulation would in general imply a fiatten- 
ing of the mass function. This is reasonable in cases where 
coagulation and collapse regulate the mass function to a 
given shape. 



APPENDIX C: 
nh-Lss PLANE 



FEEDBACK REGIMES IN THE 



The fate of the SB depends mostly on the hot phase density 
rih and the mechanical luminosity L^s. To better quantify 
the fate of the SB we restrict to the case of solar metallicity 
and /ih = 0.6, and express the hot phase density in units 
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Figure CI. Fate of SBs in the rih-Las plane. The panels indicate 
the lines separating the domains in which SBs (a) get or do not 
get into the PDS stage, (b) blow-out in the adiabatic or PDS 
stage, (c) blow out or remain confined, (d) end after exhaustion 
of SNe. 



of nii,-3 = nh/(10~ cm~ ). Moreover, H^ff is expressed in 
units of 10^ pc and the temperature of the hot phase in units 
of 10® K. Table CI reports the characteristic times of SBs. 
It is useful to analyse the fate of SBs in the nh-I/38 



plane. The condition tpds 



(ad) 



determines the possibility 



of getting into the PDS stage; it is equivalent to the line: 

(CI) 



L38 = 0.20 tIT n 



1 

h,-3 ■ 



In absence of blow-out, below this line SBs are pressure- 
confined in the adiabatic stage, above this line they can go 
into the PDS stage (Fig. Cla). The condition tpds ~ tbo' 
determines whether the blow-out is in the adiabatic or PDS 
stage; it is equivalent to: 



1-38 = 19 H, 



11/4 7/4 
off, 3 '^h,-3 



(C2) 



In case of blow-out, to the left of this line SBs blow out 
in the adiabatic stage (which is maintained even after one 
sound crossing time), to the right of it SBs blow out in the 
PDS stage (Fig. Clb). The condition tbo ~ iconf determines 
whether SBs are going to end by blow-out or confinement; 
in the adiabatic stage it is equivalent to: 



L38 = 1AH, 



S,3 



(C3) 



In the PDS stage the numerical factor is 2.9. Below this line 
SBs are pressure-confined, above it they blow out (Fig. Clc). 

Fig. C2 shows all the relations listed above and the re- 
gions they define in the nh-I/38 plane for the choice Th = 
lO'' K and -ffetr = 10'^ pc. The thick lines mark the bound- 
aries of regions where SBs end in a different way. Region (1) 
in the figure contains the SBs that blow out in the adiabatic 
stage. In region (2) SBs are confined in the adiabatic stage, 
while in region (3) SBs are confined in the PDS stage. In 



^pds 
^pda 

^conf 



2.87 X 106 l}J^^nZ^^. 



674 L 



4/11 



11 9/22^-5/11 
38 "h.-S /^h.O.e^m ^ 
7/11 1/22^-3/11 



pc 



4.12 X10« Lj/^n->^/„^.,T-ryr 



2.95 X 106 i38''^h,-3/^hA6^h,6 



5/4 



conf 



Rbo 
*bo' 

^bo' 



857 L 



1/2, 
38 ' 

1/2, 



592 1-38 n 



-1/2, 1/4 ,^-3/4 
'h,-3 f^h.O.e h,6 
1/2, 1/4 ,^-3/4 
-3 f^h.O.e h,6 



pc 
pc 



O.Oci X lU i^gg '^h,-3^h,0.6-"off,3 
r,a -1 nfi r— 1/3 l/3 l/3 rr5/3 

7.06 X lO*^ n. ' ,At, ' , yr 



r-i/3 1/3 1/3 rr5/3 
-^38 '^'■h,-3^'h,0.6-"eff,3 
10=* //,ff,3 pc 



7.06 X 10 

03 

3.11 tbo 
2.81 tbo 
1.98 flbo 
1.86 flbo 



(ad.) 
(PDS) 
(ad.) 
(PDS) 

(ad.) 
(PDS) 

(ad) 
(PDS) 
(ad) 
(PDS) 



Table CI. Characteristic times of a SB for typical values of the 
parameters. Here ^h,0.6 = /^h/O-S. 



o 




Log nh,_3 

Figure C2. As in Fig. CI, but with all lines drawn together. For 
the explanations of the five regions, see the text. 



region (4) SBs blow out in the PDS stage. Finally, in the 
closed region (5) SBs blow out in the adiabatic stage, but 
are confined before one sound crossing time. For other val- 
ues of and Th the lines move along the plane, but it is 
easy to check that the area of the triangle corresponding to 
region (5) in invariant. 

It is interesting also to consider the region of the param- 
eter space in which SBs end after all SNe have exploded. 
This happens when tbo = tufc and tconf = info. The two 
conditions give: 



-^38 
^38 



= 6.2 X 10~ 
= 1.9 X 10^ r, 



5/2 



off,3'^h,-3 
nh,-3 ■ 



(C4) 



They are both satisfied below the first Une and above the 
second (Fig. Cld). As the two Unes are parallel, this happens 
only if Heff.sT^^g^^ = 7.9. In this case a whole band in the 
plane corresponds to exhaustion. 

We conclude this discussion by noting that the plane is 
not uniformly populated, as both rih and L38 are dynamical 
variables. In other words, the probability of a certain regime 
is not determined by the area occupied in the plane but by 
the dynamics of the system. 

This paper has been typeset from a T^^X/ I^TJ^X file prepared 
by the author. 
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